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Abstract 



We consider the symmetry-breaking steady state bifurcation of a spatially-uniform equilib- 
rium solution of £'(2)-equivariant partial differential equations (PDEs). We restrict the space 
of solutions to those that are doubly-periodic with respect to a square or hexagonal lattice, 
and consider the bifurcation problem restricted to a finite-dimensional center manifold. For the 
square lattice we assume that the kernel of the linear operator, at the bifurcation point, consists 
of four complex Fourier modes, with wave vectors Ki = (a, K2 = (— /3, a), K3 = a), and 
K4 = (— a,/3), where a > /3 > are integers. For the hexagonal lattice, we assume that the 
kernel of the linear operator at the bifurcation point consists of six complex Fourier modes, also 
parameterized by an integer pair (a,/3). We derive normal forms for the bifurcation problems, 
which we use to compute the linear, orbital stability of those solution branches guaranteed to 
exist by the equivariant branching lemma. These solutions consist of rolls, squares, hexagons, 
and a countable set of rhombs, and also a countable set of planforms that are superpositions of 
all of the Fourier modes in the kernel. Since rolls and squares (hexagons) are common to all of 
the bifurcation problems posed on square (hexagonal) lattices, this framework can be used to 
determine their stability relative to a countable set of perturbations by varying a and /3. For 
the square lattice an 0(2(0; + l3) — I) truncation of the normal form is required to completely 
determine the stability of the planforms, although many of the stability results are established 
at cubic order. For the hexagonal lattice, all of the solution branches guaranteed by the equiv- 
ariant branching lemma are, generically, unstable due to the presence of a quadratic term in 
the normal form. We analyze the degenerate bifurcation problem that is obtained by setting 
the coefficient of the quadratic term to zero. We must retain terms through 0(2a — 1) in the 
normal form of the bifurcation problem. The unfolding of the degenerate bifurcation problem 
reveals a new class of secondary bifurcations on the hexagons and rhombs solution branches. 
We also analyze the bifurcation problems for E(2) + Z2-equivariant PDEs, which leads to new 
results in the hexagonal lattice case, only. 
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1 Introduction. 



Equivariant bifurcation theory [12] is a powerful tool for investigating pattern-forming instabilities 
in physical and chemical systems. This approach distinguishes between those aspects of the bifur- 
cation problem that are a consequence of symmetry and those aspects that depend on the specifics 
of the mathematical model. For example, the normal form of the bifurcation problem is derived 
using symmetry considerations alone, with details of the mathematical model appearing only in 
the numerical values of the normal form coefficients. Consequently, disparate physical and chem- 
ical systems, that nonetheless share the same symmetries, can exhibit strikingly similar behavior. 
This paper uses equivariant bifurcation theory to investigate the evolution of symmetry-breaking, 
steady state bifurcations in parameterized families of _E(2)-equivariant partial differential equations 
(PDEs), where E(2) is the Euclidean group of rotations, refiections and translations in a plane. 
Our results, which are based solely on the symmetries of the PDEs and certain features of the linear 
instability, apply to a wide variety of pattern forming systems, e.g., Rayleigh-Benard convection [3], 
models of steady cellular patterns in combustion [25] and solidification [5] and reaction- diffusion 
systems in the Turing instability regime [26]. 

Our analysis is pertinent to PDEs posed with periodic boundary conditions, and also to 
spatially-periodic states of PDEs posed on unbounded domains. The features of the linear in- 
stability that we assume are: (1) there is a time-independent, spatially-uniform trivial solution 
that loses stability in a steady state bifurcation as a parameter A is increased through the bifur- 
cation point A = Ac, and (2) at A = Ac, the trivial solution is neutrally stable to perturbations in 
the form of Fourier modes e^'^'^^''^ (k,x G R^) with finite critical wavenumber, |k| = fcc 7^ 0. A 
consequence of the rotational symmetry of the PDEs is that the kernel of the linearized problem 
at A = Ac is infinite dimensional {i.e., the length, but not the direction, of the critical wave vectors 
k is determined). One of the fundamental nonlinear problems in pattern formation is to determine 
which superpositions of Fourier modes lead to stable, steady state solutions, and to provide a simple 
characterization of these solutions, e.g. in terms of their symmetries. Our work is motivated by 
this problem. 

This paper develops a bifurcation theoretic framework for determining aspects of the stability 
of a class of spatially-periodic equilibrium solution branches that bifurcate from the trivial solution 
at A = Ac. The solutions that we investigate are periodic in two directions such that they tile 
either a square or hexagonal lattice. For example, on the square lattice, with fundamental domain 
xG[0,1)x[0,1), the steady solutions are of the form 

^^e2^«i-x + z^e"^^^^-"" + ^se'"''^'''' + z^e'^'^'-'' + c.c. + harmonics, (1.1) 

where z = (2:1,2:27-2^37-2^4) G C"*, and 

Ki = (a,/3), K2 = (-/37«)7 K3 = (/37«)7 K4 = (-«7/3)- (1-2) 

Here a > /3 > are integers, and we assume that lengths are scaled such that |Kjp = = -\- (3"^ , 
j = 1,...,4. Note that the wave vectors Ki and K2 are orthogonal, as are the wave vectors K3 
and K4, but that the angle between Ki and K3 depends on a and (3. Thus by setting two of 
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the amplitudes Zj in (1.1) to zero, we recover the bifurcation problems that address the relative 
stability of rolls and squares, or rolls and rhombs. The specific solutions that we investigate here 
are of the form 

1. Rolls (a.k.a. stripes): z = (x, 0, 0, 0), a; G R. 

2. Simple squares: z = (x,x,0,0),x G R. 

3. Two different rhombs (a.k.a. rectangles): z = (x,0,x,0) and z = (x,0,0,x), a; G R. 

4. Super squares: z = (x, x, x, x), a; G R. 

5. Anti-squares: z = (x, x, —x, —x), a; G R. 

Throughout we identify all solutions that are symmetry-related (e.g., we do not distinguish between 
rolls that are related by a rotation or translation). Note that the rolls and simple square states 
are the same, up to an overall length scaling factor, for every value of a and (3. In contrast, we 
obtain a countable set of distinct rhombs, super squares and anti-squares by varying a and (3. 
Dionne and Golubitsky [7] used the equivariant branching lemma [12, 27] to prove existence of the 
above steady solution branches for the generic _E(2)-equivariant steady state bifurcation problem. 
The equivariant branching lemma provides an algebraic criterion for existence of steady solution 
branches in steady state bifurcation problems with symmetry. Roughly, the solution branches that 
are guaranteed by the equivariant branching lemma are those that are completely characterized by 
their symmetry and by a single real amplitude; these solution branches are referred to as axial. 

Our bifurcation analysis proceeds by first restricting the space of solutions of the PDEs to those 
that are periodic with respect to some square or hexagonal lattice. Then, within this subspace 
of solutions, we invoke the center manifold theorem to reduce the bifurcation problem to a finite- 
dimensional one 

z = g(z,A), g:C^xR^C^ (1.3) 

For example, in the case of the square lattice problem (1.1) with k'^ = a'^ -\- /3^, a > /3 > 0, we 
have z G C"* and g : X R ^ C*. The vector field g inherits certain symmetries from the PDEs. 
Specifically, the group of symmetries is H-j-T^, where H characterizes the discrete symmetries of 
a fundamental domain of the lattice and is the torus of translation symmetries associated with 
doubly-periodic solutions. For the square lattice H = D4, and for the hexagonal lattice H = De. 
In the case that the PDEs are posed on a bounded domain with periodic boundary conditions, 
the periodicity of the solutions is prescribed, rather than imposed, and the symmetry is a 
consequence of the boundary conditions. 

The translation symmetry on the square lattice ensures that there are no even terms in the 
Taylor expansion of g in (1.3). However, this is not the case for the analogous bifurcation problem 
on the hexagonal lattice. The significance of this observation is that the quadratic terms in the 
Taylor expansion of g (generically) force a// of the axial solution branches to bifurcate unstably [15]. 
Our approach to investigating stable solution branches for the hexagonal lattice bifurcation problem 
is similar to that of Golubitsky, Swift, and Knobloch [13]. Specifically, we consider the following two 
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problems: (1) the degenerate bifurcation problem in which the coefficient of the quadratic term in 
the normal form is zero, and (2) the bifurcation problem for PDEs that are E(2) + Z2-equivariant, 
where the extra Z2 reflection symmetry kills the even terms in the Taylor expansion of g. Both of 
these problems arise naturally in Rayleigh-Benard convection; the degenerate bifurcation problem 
arises when the linearized operator is self- adjoint [22]; and the Z2 symmetry corresponds to a 
reflection in the midplane of the fluid layer that is present in the Boussinesq approximation [13]. 
For the hexagonal lattice, we consider solutions of the form 

^^g2^.Ki.x ^ ^^g2^.K2-x ^ ^^g2^.K3-x ^ ^^g2^.K4-x ^ ^^g2^.K5-x ^ z^g2mKe-x ^ ^ harmonics , (1.4) 

where z = (2:1,2:27-2^37-2^47-2^57-2^6) G C^. The angle between Ki and K4 is determined by an integer 
pair a > /3 > a/2 > 0. In particular, we assume that lengths have been scaled so that |Kjp = 
A;2 = q;2 _|_ ^2 _ Q,^^ j _ where 

Ki = a(0, 1) + /3(^7 -I), K4 = a(0, l) + (a - /3)(^, . (1.5) 

The wave vectors K2, K3 are obtained by rotating Ki by with K5 and Ke obtained from K4 

in the same way. In the absence of the extra Z2 symmetry, the solution branches guaranteed by 
the equivariant branching lemma are [7] 

1. Rolls: z = (x, 0, 0, 0, 0, 0), x e K. 

2. Simple hexagons: z = (x, x, x, 0, 0, 0), x £ H. 

3. Three different rhombs: z = (s, 0, 0, a;, 0, 0), z = (s, 0, 0, 0, a;, 0), and z = (s, 0, 0, 0, 0, a;), 
a; G R. 

4. Super hexagons: z = (x, x, x, x, x,x), a; G R. 

Moreover, the simple and super hexagons bifurcate transcritically, so we distinguish between the 
branch with 2:1 > and the branch with 2:1 < 0. The rolls and simple hexagons are the same state 
for every (a, (3), and there is a countable set of rhombs and super hexagons. The branching of 
super hexagons in Rayleigh-Benard convection has been investigated by Kirschgassner [16]. 

In the case that the PDEs are E(2) + Z2-equivariant, we use the equivariant branching lemma 
to show that there are five additional axial solution branches to those enumerated 1-4 above. These 
are 

1. Simple triangles: z = (ix, ix, ix,0, 0, 0), a; G R. 

2. Rhombs (called the "patchwork quilt" in [13]): z = (a;, a;, 0, 0, 0, 0), a; G R. 

3. Anti-hexagons: z = (x, x, x, —x, —x, —x), a; G R. 

4. Super triangles: z = (ix, ix, ix,ix,ix,ix), a; G R. 
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5. Anti-triangles: z = (ix, ix, ix, —ix, —ix, —ix), x £ H. 

The countable set of anti-hexagons, super triangles and anti-triangles solution branches is new. The 
simple triangles and patchwork quilt rhombs are investigated in [13]; they show that, generically, 
this branch of rhombs, composed of rolls rotated by 7r/3 relative to each other, bifurcates unstably. 

The goal of our analysis is to derive the normal form of the equivariant bifurcation problem and 
then to use it to determine the branching and linear (orbital) stability of the solutions enumerated 
above. We treat the square and hexagonal lattice bifurcation problems separately. We proceed by 
first characterizing the symmetry of each of the axial solution branches z\. These symmetries 
put restrictions on the form of the Jacobian matrix Dg evaluated on z\. Specifically, Dg(zA) 
commutes with the symmetry group Sz^- We exploit this observation to determine the eigenvalues 
of Dg and their multiplicities for each axial solution branch. We then Taylor expand the equivariant 
bifurcation problem to sufficiently high order so that the signs of the eigenvalues are determined. 
Provided certain nondegeneracy conditions are satisfied, a cubic truncation of the normal form g 
is sufficient for determining the stability of rolls, rhombs and simple squares, only. The stability 
of simple hexagons and simple triangles depend on quartic order terms in the Taylor expansion 
(or quintic order terms in the case that there is an extra Z2 symmetry). The stability of super 
squares/hexagons/triangles and anti-squares/hexagons/triangles all depend on higher order terms, 
where the order is determined by the integers a and (3 in (1.2) and (1.5). For example, we find that 
an eigenvalue (of multiplicity 2) of Dg evaluated on the super squares or anti-squares branches is 
zero unless we retain terms through 0(2(a + (3) — 1) in the Taylor expansion of g. However, we 
also find that if these two solutions are neutrally stable at cubic order, then, generically, one and 
only one of them is stable when the 0(2(a + (3) — 1) terms are taken into account. 

The bifurcation framework developed in this paper is a natural one for PDEs posed on a square 
or hexagonal domain with periodic boundary conditions. In particular, it applies when the size 
of the domain is larger than the wavelength of the instability, l/kc- For example, in the case 
of a square box of sidelength i, our bifurcation analysis applies when i is best approximated by 
\/ o? \ fP' for integers Q!,/3 G Z that satisfy a > /3 > 0. In the case that the PDEs are posed on 
an unbounded domain, our bifurcation analysis allows us to compute the stability of the periodic 
solutions in the form of rolls, simple squares, simple hexagons and simple triangles to an infinite 
number of perturbations by varying a and /3. This stability computation for the simple hexagons 
in the Benard problem is presented in [17]. 

Our paper is organized as follows. In section 2 we present a mathematical formulation of the 
bifurcation problem. We give the action of the symmetry group on the space of spatially-periodic 
solutions on the square and hexagonal lattices. We also define some terminology and outline our 
computations of the linear stability of the axial solution branches. In section 3 we characterize the 
solutions guaranteed by the equivariant branching lemma in terms of their symmetries; the neces- 
sary group theoretic computations that generate these results are banished to the appendix. The 
role of "hidden" Euclidean symmetries in the hexagonal lattice bifurcation problems is described. 
We also present some examples of the planforms associated with the axial solutions. Section 4 
contains our analysis of the square lattice bifurcation problem. We compute the eigenvalues of the 
axial solution branches in terms of the coefficients of the normal form of the bifurcation problem. 
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From this information we draw a number of conclusions about the branching and (bi)stability of 
the solutions. In section 5 we consider two bifurcation problems associated with the hexagonal lat- 
tices. We compute stability of the axial solutions for the degenerate bifurcation problem in which 
the coefficient of the quadratic term is zero. We also briefly discuss the unfolding of this bifurca- 
tion problem and present an example bifurcation diagram that indicates the secondary bifurcation 
points on the axial solution branches. We then consider the bifurcation problem in the case that 
there is an extra Z2 symmetry. Section 6 contains our conclusions. 

2 Problem Formulation. 

We consider parameterized families of partial differential equations which we write in evolutionary 
form, 

^u(x,0 = F(u(x,0,A), (2.1) 

where FiA'xR^J^isa nonlinear operator between suitably chosen function spaces, X and y, 
and A G R is the bifurcation parameter. Here u : X R ^ R" is a function in A' of a spatial 
variable x G R^ and time t. For simplicity we have suppressed any possible dependence of u on 
a third bounded spatial variable y. The y-dependence is important to our analysis only in so far 
as it can introduce additional symmetry to the problem. In this paper, we consider only the case 
where the additional symmetry is a reflection. 

2.1 Symmetries of the PDEs. 

We assume that (2.1) has Euclidean symmetry. The Euclidean group E{2) is the group of motions 
m R2 that preserve distances, i.e. rotations, reflections and translations. We denote elements of 
E{2) by (/i, d) where h G 0(2) is an orthogonal transformation (a reflection or rotation) and d G R^ 
is a translation. The action of (/i, d) on x G R^ is defined by 

(/i, d)x = /ix d . (2.2) 

This action of E{2) on R^ forces the product of (/ii,di) and {h2.,d2) to be defined by 

(/ii, di)(/i2, d2) = (/ii/i2, di + hrd2) . (2.3) 

Hence E{2) is the semi-direct product (denoted by -(-) of the groups of orthogonal transformations 
and translations; specifically, E{2) = 0(2)-i-R^, where R^ is a normal subgroup of E(2). 
The Euclidean group acts on the vector- valued function u : R^ X R ^ R" as follows 

7.u(x, t) = A^u(7-^x, t) (2.4) 

for all 7 G E(2), where is an orthogonal n X n matrix. For example, when (2.1) is a system 
of n reaction- diffusion equations, = I„ for all 7 G E(2), where I„ is the n X n identity matrix. 
For the Navier-Stokes equations in the plane, A^ = I2 when 7 is a translation, and A^ is a 2 X 2 
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orthogonal matrix when 7 = (h,0), h G 0(2). Our assumption that (2.1) has Euclidean symmetry 
means that F is _E(2)-equivariant, i.e., 

7.F(u(x,0,A) = F(7.u(x,0,A) (2.5) 

for all 7 G E(2), where the action of 7 on the vector- valued function F is given by (2.4). 

The symmetry of the problem is enlarged from E(2) to E(2) + Z2 for some of the motivating 
applications. For example, in certain Rayleigh-Benard convection problems Z2 is a reflection in the 
mid-plane of the fluid layer [13]. Reaction- diffusion systems, 

^ = PV2u + g(u), g(-u) = -g(u), (2.6) 
where P is a (constant) matrix of diffusion constants, also possess a reflection symmetry u — u. 

2.2 Linear analysis and the symmetry-breaking bifurcation. 

We assume that there is a Euclidean-invariant time-independent solution of (2.1) for all values of 
A. This corresponds to a spatially uniform equilibrium, which, without loss of generality, we take 
to be u = 0. We assume that this trivial solution is linearly stable for A < Ac, and linearly unstable 
for A > Ac. Moreover, we assume that there is a symmetry-breaking steady state bifurcation at 
A = Ac. At this bifurcation point the zero solution is neutrally stable to perturbations in the form 
of spatial Fourier modes e^'^'^^''^ with k G R^, |k| = kc- The neutral stability curve, A = A(|k|), is 
determined by seeking equilibrium solutions u = uj^e^'^*'^'-''- of the linearization ^ = L^u of (2.1) 
at u = 0. (Here U]j is a constant ra-dimensional vector.) We refer to the equilibrium solutions 
^j^g27r8k-x^ |k| = kc, of the linearized problem at A = Ac as the critical or neutral modes. A typical 
neutral stability curve is depicted in Figure la. Note that the minimum of the neutral curve 
occurs at (|k|,A) = (kc,Xc), where kc is nonzero and finite. A consequence of the 0(2) C E(2) 
symmetry is that the minimum corresponds to a circle of radius kc in the two-dimensional k-space 
(see Figure lb). Without loss of generality, we set Ac = in the remainder of the paper. 

2.3 Spatially doubly-periodic solutions. 

We restrict our bifurcation analysis to solutions u(x, t) of (2.1) that are doubly-periodic with respect 
to some square or hexagonal lattice C Specifically, the planar lattice C is generated by two linearly 
independent vectors ii,i2 G R^, i.e., 

£ = {m^i +n24 G R^ : rai,ra2 G Z}. (2.7) 

In this paper we consider two cases that satisfy \ii\ = \i2\. (1) the square lattice with 

it = (1,0), 4 = (0,1), (2.8) 

and (2) the hexagonal lattice with 
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We say that a function u(x, t) is C — periodic if 



u(x + ^,i) = u(x,i) for all i e jC. (2.10) 
Moreover, we assume that ^-periodic solutions of (2.1) can be expressed in a Fourier series 

Uj{^,t)= (^i,- k(0 e'"*-^ + c.c.) , j = l,...,n, (2.11) 

where Uj i^ G C is the time-dependent amplitude of the k*'' Fourier mode. The wave vectors k lie in 
the dual lattice to £, denoted C* . Specifically, C* is generated by two linearly independent vectors 
ki,k2 G R^, where k^ ■ ij = Sij (the Kronecker delta): 

jC* = {niki + n2k2 G : ni, G Z}. (2.12) 

An important consequence of restricting the solution space of (2.1) to ^-periodic functions is 
that the spectrum of the linear operator La is rendered discrete. Hence, we expect the center 
manifold theorem [19] to apply at the bifurcation point. Specifically, for the problems of interest, 
this restriction ensures that there are only a finite number of zero eigenvalues at the bifurcation 
point, with all other eigenvalues bounded away from the imaginary axis. The dimension of the 
bifurcation problem depends on the number of points k G that lie on the critical circle of radius 
kc- For the square lattice, we consider the case where the critical circle intersects 8 points in £*, 
and, for the hexagonal lattice, we consider the case where the center manifold is 12- dimensional (see 
Figure 2). We note that the simpler cases where 4 points intersect the critical circle for the square 
lattice, and 6 points for the hexagonal lattice have already been analyzed (see [12], and references 
therein); we recover many of the results of these earlier studies in the course of our analysis. 

In what follows we identify the kernel of the linear operator Lq, 

s 

ker(Lo) = {u = ^ ^je^'^^'^^-^Uj + c.c. : Zj G C, \Kj\ = k^} , (2.13) 
with the vector space 

s 

V = {v = J2 Zje^^'^'-'' + c.c. : z^ G C, |Kj| = A;J = . (2.14) 

In (2.13) Uj is a constant ra-dimensional vector associated with the Kj Fourier mode. The isomor- 
phism between V and is defined by 

V ^ z = izi,zi, . . .,Zs) . (2.15) 

As a vector space over the reals, dim(F) = 2s. As mentioned above, this paper focuses on the case 
s = 4 for the square lattice and s = 6 for the hexagonal lattice. 
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The PDEs, restricted to the center manifold, lead to a system of ordinary differential equations 

z = g(z,A), g:C"xR^C^ (2.f6) 

Here ^(0, A) = and the Jacobian matrix at the bifurcation point, Dg(0, 0), is the zero matrix. In 
the next section we describe the symmetries inherited by the bifurcation problem from the PDEs. 
In particular, if T is the symmetry group of the bifurcation problem (2.16), then g(z. A) satisfies 
the usual equivariance condition 

7g(z,A) = g(7z,A), forall7Gr. (2.17) 
2.4 Symmetry of the restricted bifurcation problem. 

The symmetry of the PDEs (2.1), reformulated in the space Xc of ^-periodic functions, is a compact 
group r. Specifically, T is the largest group, constructed from E(2), that preserves Xc, i.e., 
j.Xc C Xc for all 7 G r. As with E(2), T has a semi-direct product structure, namely T = H-j-T^, 
where H C 0(2) is the finite group of rotations and refiections that preserve the lattice and 

~ 'R? /C is the torus of translations. The discrete group H is called the holohedry of the lattice; 
in the case of the square lattice, H = D4, while H = De for the hexagonal lattice. (Recall that D„, 
the dihedral group of order 2ra, is the group of symmetries of a regular ra-gon.) In this paper we also 
consider the case where T is enlarged to P + Z2. In the remainder of the paper, let P^ = D4-i-T^ 
and Th = De+T^, while P, without a subscript, refers to Ps(-|-Z2) and/or Thi + Z2)- 

We note that if the PDEs (2.1) are posed with periodic boundary conditions in two linearly 
independent directions, then the symmetry of the full problem is exactly P = H-j-T^. In this situ- 
ation, the symmetry is a consequence of the periodic boundary conditions, and H characterizes 
the symmetries of the spatial domain. 

In this paper, we assume that the action (2.4) of 7 G P on vector- valued functions u(x,i) leads 
to an action of P on real- valued scalar functions v(x,t) in (2.14) given by 

7.^;(x,0 = ^^(7"'x,0 • (2.18) 

This is the case for all of the applications mentioned in the introduction. (See [1] for examples of 
"pseudo-scalar" PDEs that do not meet this criterion.) 

Square lattice case. 

For doubly-periodic solutions on a square lattice we take the generators of the dual lattice C* to be 

ki = (1,0) and = (0,1) . (2.19) 

Thus the wave vectors k G in (2.11) have the form (rai,ra2), where rii and n2 are integers. 
Moreover, we assume that lengths in the original PDEs have been scaled so that kc = \/ o? \ fP' 7^ 
for some integers a and /3. Alternatively, we could hold fixed and scale the lattice C 
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Table 1: Translation Free (absolutely) Irreducible Representations For the Square Lattice. 



dim(F) 


K's 


4 

(s=2) 


Ki = ki = (l,0) 
K2 = k2 = (0,1) 


8 

(s=4) 


Ki = aki + /3k2 = (a,/3) 

K2 = -/3ki + ak2 = (-13, a) 

K3 = /3ki + ak2 = ip,a) 

K4 = — aki + /3k2 = (—a, (3) 

a, (3 e Z, a > (3 > 0, (a, (3) = 1 *, 

a and (3 are not both odd. 



* (a, (3) = 1 means that a and (3 are relatively prime. 



The relevant representation of the symmetry group F^ = D4-i-T^ is determined by considering 
its action on the complex amplitudes Zj of the critical Fourier modes in (2.14). The irreducible 
representations of F^ are either 4-dimensional or 8-dimensional, in which case there are two or four 
complex Fourier amplitudes, respectively (i.e., s = 2 or s = 4). Examples of these two different 
cases are depicted in Figure 2a for kc = I and kc = i.e., for {ct,/3) = (1,0) and {ct,/3) = (2, 1). 
Note that it is also possible for the critical circle to intersect more than eight points in the dual 
lattice, e.g., [{ kc = 5 then there are four (real) Fourier modes associated with (a, (3) = (5,0) and 
eight associated with (a, (3) = (4,3). We do not consider these special cases here. (See Crawford 
[6] for an application of these higher- dimensional reducible representations.) 

Following Dionne and Golubitsky [7] we require the representation of F^ to be not only irre- 
ducible, but also translation free. A representation is translation free if there are no (non-trivial) 
translations in F^ that act trivially on (2.14). This requirement ensures that we have found the 
finest lattice C that supports the neutral modes (2.14) [7]. Table 1 gives the values of the criti- 
cal wave vectors for the translation free (absolutely) irreducible representations, henceforth simply 
called representations. Note that there is just one four- dimensional representation. It is the one 
that applies when the periodicity of functions in Xc coincides with the wavelength of the instability, 
i.e. kc = |ki| = |k2|. The bifurcation problem associated with this representation of F^ has been 
studied extensively. In this paper we focus on the eight-dimensional representations associated with 
the integer pairs (a, (3) where a > (3 > (see Figure 3). The additional requirements in Table 1, 
namely that a and (3 are relatively prime and not both odd, ensure that the representation is 
translation free, and hence that the set of all critical modes (2.14) cannot be supported by a finer 
lattice C (see [7]). 

D4 C Ts is generated by a counterclockwise rotation i?^/2 by 7r/2 about the origin and a 
refiection t^^ through the si-axis. The elements of C F^ are denoted by = {61,62), where 
^1, ^2 G [0, 1). The action of F^ on V given by (2.14) with s = 2 in Table 1 induces an action of F^ 
on generated by 

i2,/2(z) = (^2,^i), (2.20) 
r,,(z) = (^1,^2) (2.21) 
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and 



0(z) = (e-2-ki-0^i,e-2-k2-e^^^ 

= (e-^^'^'zi,e-^^'^^Z2) . (2.22) 

The action of on V for s = 4 in Table 1 induces an action of on generated by (cf. Figure 3) 

R7t/2{z) = {Z2,ZI,Z4,Z3) , (2.23) 

r^i(z) = (^4,^3,^2,^1) , (2.24) 

and 

0(z) = (e-2^''Ki-©^i,e-2^''K2-©^2,e-2^''^3-®^3,e-2^''^*-®^4) (2.25) 
Hexagonal lattice case. 

For doubly-periodic solutions on a hexagonal lattice the generators of the dual lattice C* are 

ki = (0,l) and = (^3/2, -1/2) . (2.26) 

We assume that lengths in the original PDEs have been scaled so that kc = \/ o? + f]"^ — a(3 ^ 
for some integers a and [i. 

The relevant representation of the symmetry group = De+T^ is determined by considering 
its action on the complex amplitudes of the critical Fourier modes at the bifurcation point. As for 
the square case, the neutrally stable modes at A = are given by (2.14). In this case the irreducible 
representations of are either 6-dimensional or 12-dimensional, in which case there are three or 
six complex amplitudes, respectively {i.e., s = 3 or s = 6). Examples of these two different cases 
are depicted in Figure 2b for kc = 1 and kc = V^, i.e., for {a,f3) = (1,0) and {a,f3) = (3,2). The 
values of the critical wave vectors for the translation free (absolutely) irreducible representations 
are summarized in Table 2. Note that there is just one six- dimensional representation which is 
associated with the case where the periodicity of functions in Xc coincides with the wavelength of 
the instability, i.e. kc = |ki| = |k2|. The bifurcation problem associated with this representation 
of Th has been studied extensively [4, 13]. In this paper we focus on the twelve- dimensional 
representations associated with the integer pairs (a, (3), a > (3 > a/2 > (see Figure 4). The 
restriction (a, (3) = (3,a + (3) = 1 ensures that the representations are translation free. 

De C Th is generated by a counterclockwise rotation by 7r/3 about the origin and a 

reflection t^^ through the si-axis. The elements of C F/^ are denoted by = ^li^i -|- ^2^2, where 
= (I/V3, 1), £2 = (2/V3,0), and 6*1, 6*2 G [0,1). The action of Th on V given by (2.14) with 
s = 3 in Table 2 induces an action of Th on C"^ generated by 

i2,/3(z) = (^3,^1,^2), (2.27) 
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Table 2: Translation Free (absolutely) Irreducible Representations For the Hexagonal Lattice. 



dim(F) 


K's 


6 


Ki = ki = (0,l) 


s=3 


K2 = k2 = (73/2,-1/2) 




Kg = -ki -k2 = (-73/2,-1/2) 


12 


Ki = aki + /3k2 


s=6 


K2 = {-a + /3)ki - Q!k2 




K3 = -/3ki + {a- /3)k2 




K4 = aki + (a — /3)k2 




K5 = -/3ki - Q!k2 




K6 = (-a + /3)ki + /3k2 




a,l3 e Z, a > l3 > a/2 > 0, 




(q!,/3) = 1 and {3,a + f3) = l *. 



* (3, a + /3) = 1 means that a + /3 is not a multiple of 3. 



r,,(z) = (^1,^3,^2) (2.28) 

and 

0(z) = (e-2-ki-0^i,e-2-k2-e^^^g-2.k3.e^^^ 

= (e-^^'^^zi,e-^^'^'z2,e^^'^^^+^'h3) . (2.29) 
The action of Th on F for s = 6 in Table 2 induces an action of Th on generated by (c/. Figure 4) 

R^/siz) = (^2,^3,^1,^5,^6,^4) , (2.30) 

r^l(z) = (z6,Z5,Z4,Z3,Z2,Zi) (2.31) 



and 



0(z) = (e-2-Ki-0^i, e-2-K2-e^^^ g-2..K3.e^^^ g-2..K4.e^^^ g-2..K5.e^^^ g-2..Ke.e^^^ 

p^-27Ti{a0t+{a-fj)02) p^-27Ti{-fj0t-a02) p^-27Ti{{-a+fj)0t+fj02) ^^-^ _ (2.32) 



Additional Z2 symmetry. 

In this paper we consider the possibility that there is an additional Z2 symmetry so that the 
bifurcation problems are equivariant with respect to (H-j-T^) + Z2, where H = D4 or H = Y)q. 
We assume that k G Z2 takes v to —v, where v £ V is given by (2.14). This induces the following 
action on z G C*: 

k(z) = -z . (2.33) 
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The additional reflection symmetry has no effect on the bifurcation problems associated with 
the square lattice. This observation, for the four- dimensional representation of (D4-i-T^) + Z2, is 
made in [24]. The case of the eight-dimensional representations in Table 1 is the same. Speciflcally, 
we note that the translation (^,^) G in (2.25) acts on z in the same way as the reflection k 
in (2.33). Hence, we need consider the effect of the additional reflection for bifurcation problems 
associated with hexagonal lattices only. 

2.5 Overview of calculations. 

The goal of this work is to compute stability of particular solution branches of the bifurcation 
problem (2.16). In particular, we compute the stability, at bifurcation, of those solutions that are 
guaranteed to exist by the equivariant branching lemma [12]. This lemma provides an algebraic 
criterion for existence of solution branches associated with particular subgroups of T. Speciflcally, 
we specify the symmetry of an equilibrium solution z G by the isotropy subgroup Sz C T, where 

Sz = {(7 G r : dz = z} . (2.34) 

A subgroup S C r is an isotropy subgroup if there exists a z G for which Sz = S. Associated 
with each isotropy subgroup S C T is a vector subspace of C*, called the flxed point subspace and 
denoted Fix(S), where 

Fix(S) = {z G : dz = z, for all d G S} . (2.35) 

The equivariant branching lemma states that provided certain (generic) conditions are satisfled 
by the bifurcation, there exists a branch of equilibrium solutions, bifurcating from the origin at 
A = 0, with symmetry S for each isotropy subgroup S C T that satisfles dim(Fix(S)) = l (see [12]). 
In the next section we determine, up to conjugacy, all isotropy subgroups for the 8-dimensional 
representations of D4-i-T^ and the 12- dimensional representations of De-j-T^ and (De-j-T^) -|- Z2. 
We then depict examples of those planforms guaranteed to exist by the equivariant branching 
lemma. Following [11], we refer to isotropy subgroups with 1-dimensional flxed point spaces as 
axial and the associated spatially doubly-periodic solutions as axial planforms. 

To determine the linear stability of an equilibrium solution branch z\ of the F-equivariant 
bifurcation problem z = g(z. A), we compute the eigenvalues of the Jacobian matrix Dg(zA, A). A 
simple consequence of the equivariance (2.17) of g : X R ^ is that the Jacobian matrix 
evaluated at z = za commutes with all a G ^zx- This condition puts restrictions on the form of 
Dg(zA,A). We exploit this observation to determine the eigenvalues of the Jacobian matrix for 
each axial subgroup. We then determine which of the eigenvalues are forced, by the continuous 
translation symmetry T^, to be zero. Finally we determine expressions for the nonzero eigenvalues 
in terms of the leading order terms in the Taylor expansion, at z = 0, of the general smooth 
F-equivariant vector held g(z,A). These results determine conditions, for stability of the axial 
planforms, that must be satisfled by the coefficients of certain terms in the Taylor expansion of g. 

In the case of the hexagonal lattice with De-j-T^ symmetry, the Taylor expansion of the general 
equivariant vector field possesses a term that is quadratic in the Zj. It then follows, from a theo- 
rem due to Ihrig and Golubitsky [15], that generically the axial planforms all bifurcate unstably. 
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Therefore, in order to investigate stable solutions within the setting of a local bifurcation problem, 
we study the degenerate bifurcation problem obtained by setting the coefficient of the quadratic 
term to zero. We also briefly indicate some of the interesting possibilities for secondary bifurcations 
within an unfolding of this degenerate bifurcation problem. 

3 Group Theoretic Results. 

3.1 Isotropy subgroups. 

In this paper we follow the convention of identifying all solution branches that are on the group 
orbit Tz\ of a particular branch z\. Thus we classify isotropy subgroups of T by conjugacy class 
since the isotropy of a point z\ G is conjugate to the isotropy of a point on its group orbit. 
Specifically, S-yz^ = T^z^T"^- (Recall that two subgroups Si, S2 C T are conjugate if S2 = 7X17"^ 
for some 7 G T.) 

In this section, we determine, up to conjugacy, all isotropy subgroups and their fixed point 
spaces for the 8-dimensional representations of F^, and for the 12- dimensional representations of Th 
and F/i + Z2. In each case we summarize the results in a table; the computations that are necessary 
to generate the tables are relegated to the appendix. 

Square lattice case. 

We list in Table 3 the isotropy subgroups, up to conjugacy, of acting on together with their 
generators, their fixed point subspaces and the dimensions of the fixed point subspaces. Note that 
the pure translation subgroups, denoted S^, 81^27 ^1,3 and 81^4, depend on the values a and (3 and 
hence are not the same for all 8-dimensional representations. Associated with these fixed point 
subspaces are planforms in V that are periodic with respect to a finer lattice than C There are 
six isotropy subgroups for which dim(Fix(S)) = l; these are the axial subgroups. The equivariant 
branching lemma guarantees the (generic) existence of a branch of equilibria associated with each 
of these subgroups S. 

Hexagonal lattice cases. 

Table 4 contains the isotropy subgroups, up to conjugacy, of Th acting on together with their 
generators, their fixed point subspaces and the dimensions of the fixed point subspaces. The isotropy 
subgroups associated with Th + Z2 are listed in Tables 4 and 6. The equivariant branching lemma 
applies to six of the isotropy subgroups in the case that F = F/^ and ten in the case that there is an 
extra Z2 symmetry. Note that some of the isotropy subgroups in Table 4 are modified by the extra 
Z2 symmetry as indicated in Table 5. In this table we denote elements of Th = De-j-T^ by (h, 0) 
and elements of Th + Z2 by ((h,Q),Id) and ((/i,0),k). Here h G De, G and Z2 = {Id,K}, 
where Id specifies the identity element of a group. 
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Table 3: Isotropy Subgroups S (up to conjugacy) of F^. 







Fix(S) 


T)im(Fix(T,]) 


r 


B In T 


z = (2:1,2:27^^37^^4) = 







n /o 

^^TY 1 2^ ' Xl 




Zl = Z2 = Z3 = Z4 e 


1 




-"-71/27 'v'a;i7 ^2' 2// 




Zl = Z2 = -Z3 = -Z4 e R 


1 


Do 


-"■717 'a 




Zl = Zs, Z2 = Z4 E 


2 


D5 
^2 


/? r 

-^t'7r7 ' Xl 




Zl = Z4, Z2 = Zs E 


2 


fix 


-"-717 V 7 V 2 ' // 




= —Z4, Z2 = Z3 £ H if a is odd. 


2 






Zl 


= Z4, Z2 = —Z3 G R if /3 is odd. 


2 


Df 


-^717 {Txi 7(272)) 




2^1 = — 2:47 ^2 = ~^3 G R 


2 


Z4 


-^71/2 




2^1 = 2^2 7 ^3 = ^4 G R 


2 


zl 






^1 = ^3 7 ^2 = ^4 G C 


4 


Zf 






Zl = Z4, Z2 = Z3 E C 


4 


Zf 


(^.17(^0)) 


Zl 


= —^4, 2^2 = ^3 G C if a is odd. 


4 






Zl 


= ^4, 2^2 = —^3 G C if /3 is odd. 


4 


Z2 






^l7 ^27 ^37 2:4 G R 


4 


1 


Id 




^l7 ^27 ^37 Z4 £ C 


8 


z^+si 


R.. si 


Zl e B., Z2 = Zs = Z4 = 


1 


Z4 + Si^2 


-^71/27 Si, 2 ^''^ 




Zl = Z2 e B., Zs = Z4 = 


1 


Z2 + Sl,2 


-^717 Si, 2 




zi,Z2 eB., Z3 = Z4 = 


2 


D^+Si,3 


i?7r7 Si,3 




Zl = Zs e B., Z2 = Z4 = 


1 


^2 + ^1,3 


R-K^ Si,3 




Zi,Z3 eB., Z2 = Z4 = 


2 


Df+Si,4 


r^i, Si,4 (-^^ 




Zl = Z4 e B., Z2 = Zs = 


1 


+ 


i?^, Si,4 




Zi,Z4 eB., Z2 = Z3 = 


2 



the generators of are given in equations 2.30-2.32. 
^) i?^ = R'^12'1 Td = '''xiR'l/2 ^ reflection through the line containing ii + ^^2 (c/. Fig. 3). 

= {(/3s, -as) G T2 : s G R} . 
^) Si,2 is generated by (^, ^), (^, ^) G T^. 
^) Si,3 is generated by (^, ^), (^, ^) G T^. 
^) Si,4 is generated by (2^, Jp) ^ T\ 
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Table 4: Isotropy Subgroups S (up to conjugacy) of Th and F/^ + Z2. 





Generators of S 


Fix(S) 


Dim(Fix(i;)) 




R /, (k) 


Z = r ^1 Zo Z'x Za Za \ — ^ 





Dr 


i? /Q r 


P'l — P'o — P'q — P'^ — P'c — Za ^ 'W. 


1 




-^7r/3' 


Z-\ — P'o — P'q — ~Z A — "Fc — ~Z a C^j 


2 


-'--'3 


-^7r/3' '^1 


^1 — ^2 — ^3 — ^4 — ^5 — ^6 t ^ 


2 




7? , 

-ft7r/3 


^1 — ^2 — ^3? ^4 — ^5 — ^6 t ix 


9 
z 




7? T- 


"y.. — "y^ "y^ — "y— "y^^ — "y . ^ l-< 

^1 — ^67 ^2 — ^5 7 ^3 — ^4 t ix 


Q 


Z3 




^1 = ^2 = ^3, ^4 = ^5 = ^6 G C 


4 


z^ 




Zi = 2:4, 2:2 = ^67 2^3 = ^5 G C 


6 


Zf 




^1 = ^67 ^2 = ^5, ^3 = ^4 G C 


6 






^17^27^37^47^57^6 G R 


6 


1 


Id 


^17^27^37^47^57^6 G c 


12 


z^+si 


R.. Si 


2^1 G R, 2^2 = ^3 = ^4 = ^5 = ^6 = 


1 


D-+Si,4 


Rtvi T'ai Si, 4 


2:4 = 2:4 G R, 2:2 = 2:3 = 2:5 = 2:5 = 


1 


Z2 + Sl,4 


-Rtt, Si,4 


2:4 , 2:4 G R, 2:2 = 2:3 = 2:5 = 2:6 = 


2 




p _ (d) c 
-^7r? 'm ? "^1,5 


2^1 = 2:5 G R, 2:2 = 2:3 = 2:4 = 2:5 = 


1 


Z2 + Sl,5 


-Rtt, Si, 5 


2^1 , 2:5 G R, 2:2 = 2:3 = 2:4 = 2:6 = 


2 


Df+Si,6 


Rti") ^xi •} Si, 6 


zi = zq a R, 2:2 = 2:3 = 2:4 = 2:5 = 


1 


Z2 + Sl,6 


-Rtt, Si,6 


2^1 , 2:6 G R, 2:2 = 2:3 = 2:4 = 2:5 = 


2 


Z6 + Si^2,3 


-^71/37 Si, 2,3 


2^1 = 2:2 = 2:3 G R, 2:4 = 2:5 = 2:5 = 


1 


Z3 + Si^2,3 


^^/3' ^1,2,3 


2:4 = 2:2 = 2:3 G C, 2:4 = 2:5 = 2:6 = 


2 


Z2 + Sl,2,3 


-Rtt, Si, 2,3 


^l7 ^2 7 ^3 G R7 ^4 = ^5 = ^6 = 


3 


Si, 2,3 


Si, 2,3 


^1 7 ^2 7 ^3 G C , 2:4 = 2:5 = 2:6 = 


6 



the generators of are given in equations 2.23-2.25, and k G Z2 is given in (2.33). 
Tn = R^i^Tx^ is a reflection through the line containing the vector li — 2I2 (cf. Fig. 4). 
The generators of S^, Si, 4, Si, 5, Si,6, and Si, 2, 3 are given in Table 5. 
(''^ Tm = R^i'iTxi is a reflection through the line containing the vector li + 12 {cf. Fig. 4). 
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Table 5: Generators of Subgroups and S* of Th and F/^ + Z2. 



sVs* 


Generators of S^, S* C Th 


Generators of S^, S* C T/^ + Z2 


si 


{ 1 ri iH Q f -1 /'V e Tt-ii TxrhoTO e rt 

[1 iJi^ j-J-o \.\ iJL-o M-2 ) •) W lit;! t; -o xz -l-V 


Uj ^ jj -CI ^ ^ 2^ (-Xo J -1,2 / ? W lit:! t: o t XV 


Sl,4 


(Id, j:L, ^1 «-2, ^2), 

(Id, a^-U ^1 + «-2/3 ^2) 


((Id, 2(«2-2a/3) ^i)' '^)' 

((^^' 27^ (^1 - 2^2)), 


Si, 5 


(Id, h ^2), 

(Id, ^2_p2 f-l „2_^2 


((^^' 2(^(^1-^2)), 
aid, ^(^1+4)), k) 


Sl,6 


(Id, ii + 2„7_% ^2), 

(Id, 2a-fj 1 2af3-f3^ ^2) 


((^^'2(2^(2^1-^2)), k). 


Si, 2,3 






(-^^' a2-a/3+/32 „2-a/3+/32 ^2), 


((^^' „2-„^ + ,2 h „2_"J,V ^2), ^^), 


(Id, a^-afj+fj^ „2_„^_,_^2 -^2) 


({Id, „2_„^ + ^2 „2_„^ + ^2 ^2), Id) 


Si, 2 


not applicable 


((Id, 2(a^-af3+f3^) 2(a^-af3+f3^) ^2)7 «), 
((-^^' 2(a2-a/3+/32) 1 2{a^-af3+f3^) ^2)7 '«) 



Table 6: Isotropy Subgroups S (up to conjugacy) of Th + Z2. Also see Table 4. 



s 


Generators of S 


Fix(S) 


Dim(Fix(S)) 


De 


((^7r/3,0),K),((r^l,0),K) 


-^1 = -^2 = -^3 = --^4 = -Z5 = -Zq e Ri 


1 


De 


((i2,/3,0),K),((r,,,0),/rf) 


-^1 = -^2 = -^3 = -^4 = -^5 = -^6 G 


1 


De 


((i2,/3,0),/rf),((r,,,0),K) 
((^'/3'0),/rf),((r.„0),K) 


^1 = ^2 = ^3 = ~^4 = = G R 


1 


Ds 


^1 = ^2 = ^3 = —^4 = —^5 = —^6 G C 


2 


D3 


((^'/3'0),/rf),((r.,0),K) 


^1 = ^2 = ^3 = —^4 = —^5 = —^6 G C 


2 


Ze 


((i2,/3,0),K) 


^1 = ^2 = ^3? ^4 = ^5 = ^6 G Ri 


2 


D2 


((t,,,0),k),((R^,0),k) 


^1 = —^67 ^2 = —^57 ^3 = —^4 G Ri 


3 


D2 


((T,„0),K),((R^,0),Id) 


^1 = —Ze, Z2 = —Z5, Z3 = —Z4 G R 


3 


D2 


((T,„0),Id),((R^,0),K) 


^1 = Ze, Z2 = Z5, Z3 = Z4 E Hi 


3 


Z2 


((r^i,0),K) 


Z\ = —Zq, Z2 = —Z5, Z3 = Z4 E C 


6 


Z2 


((r„,0),K) 


Z\ = —Z4, Z2 = —Zq, Z3 = Z5 E C 


6 


Z2 


((R^,0),k) 


z\, Z2, z^, 2:4, 2:5, zq G Ri 


6 


Z6 + Sl,2,3 


((i2,/3,0),K),Si,2,3 


-^1 = -^2 = -^3 G Ri, 2:4 = 2:5 = 2:6 = 


1 


Z2 + Sl,2,3 


((_R^,0),k), Si,2,3 


Zl,Z2, Z-i G Ri, 2:4 = 2:5 = 2:6 = 


3 


Z2 + Sl,2 


((i2^,0),/rf),Si,2 


^1 , ^2 G R, ^3 = ^4 = ^5 = ^6 = 


2 



(a) ^j^g generators of Si, 2 and Si, 2, 3 are given in Table 5. 
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3.2 Hidden symmetries. 



In the case of the hexagonal lattice, we must take certain hidden symmetries into account in our 
analysis. The hidden symmetries are elements of the full Euclidean group that are not present when 
the PDE is restricted to a periodic domain. Nonetheless, these hidden symmetries are manifest 
in certain fixed point subspaces of and put restrictions on the normal form of the bifurcation 
problem (2.16). We refer the reader to Crawford [6] for a detailed treatment of hidden Euclidean 
symmetries in mode interaction problems. 

For the 12- dimensional representation of Th, the hidden Euclidean symmetries are generated 
by a refiection through the line containing the vector (3ii — ai2, denoted by t^^ G E{2). This 
refiection acts on z GFix(Si^2,3) — {(-2^17-2^27-2^37 0,0,0) : Zi G C} in the same manner that t^^ G De 
acts on the six- dimensional representation of (cf. equation 2.28). Specifically, 



Once the hidden refiection t^^ G E(2) is included we can reformulate the F/i-equivariant bifurcation 
problem (2.16) restricted to the subspace {(2:1, 2:27 -2^87 0, 0, 0) : Zi G C} as a bifurcation problem for 
the six- dimensional representation of Th in Table 2. This is consistent with the observation that 
when the three amplitudes Z4,Z5,ze are zero, the solutions are all periodic with respect to a finer 
lattice; the basis vectors for the dual to this finer lattice are Ki and K2. 

Only in the case of Th + Z2 does the inclusion of the hidden refiection ensure the existence of an 
additional axial planform. Specifically, the line containing z = (1, 1, 0, 0, 0, 0) in the two-dimensional 
fixed point space Fix(Z2-(-Si^2) — {(2^1, a;2, 0, 0, 0, 0) : Xi G R} is fixed by TxiR-k/?,- 

For the 8-dimensional representation of F^, there is a hidden refiection that acts on z G 
{(2:1, 2:27 0, 0) : Zi G C} as follows: 



However, including this hidden symmetry does not change the normal form of the bifurcation prob- 
lem nor does it lead to new axial planforms. (Note that every point in the subspace {(2:1, 2:2, 0, 0) : 
Zi G C} is on the group orbit of one in Fix(Z4-i-Si^2) — {(2^1, a;2, 0, 0) : Xi G R}, and that the hidden 
refiection t^i acts trivially on points in Fix(Z4-i-Si^2)-) 

3.3 Axial planforms. 

In this section we present examples of the planforms associated with the axial subgroups. Specifi- 
cally, we present grey scale plots of the function f(x) in (2.14) for a representative point z GFix(S). 
We do this for each conjugacy class of subgroups S C F that fix one- dimensional subspaces. 

In the case of the square lattice, there are six conjugacy classes of isotropy subgroups that fix 
one- dimensional subspaces of C*. These are listed in Table 7. Figure 5 presents the associated 
axial planforms v in the case that {a,[i) = (2, 1). The rhombic, super square and anti-square states 
depend on (q!,/3). The rolls and simple squares are, up to scaling of the spatial variable x, the same 
for each 8-dimensional representation of F^. 



T"2;i(^^l7^^27^^37 070,0) = (2:1, 2:3, 2:2, 0,0,0). 



(3.1) 



T"2;i(^^l7^^27 07 0) = (2:1,2:2,0,0). 



(3.2) 
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Table 7: Axial Planforms for 8- dimensional representations of F^. 



Name 


Nomenclature 


representative point z G Fix(S) C C"* 


Rolls 


R 


z = (1,0,0,0) 


Simple Squares 


SiS 


z = (l, 1,0,0) 


Rhombs (^) 


Rhsl,a,/3 


z = (l, 0,1,0) 


Rhombs 


Rhs2,a,/3 


z = (l, 0,0,1) 


Super Squares 


SuS„,/3 


z = (l, 1,1,1) 


Anti-squares 


AS„,/3 


z = (l, 1,-1,-1) 



(«) This state is also called "stripes". 
(^) Rhombs are also called "rectangles". 

(^•^ Super squares are called simply "squares" in Dionne and Golubitsky [7]. 



In the hexagonal lattice case there are six different axial planforms for the twelve- dimensional 
representations of F = F/^ and eleven for F = F/^ -|- Z2, including the state Rh/io determined by 
taking into account the hidden reflection symmetry. These are listed in Table 8 and depicted in 
Figures 6 and 7. In the case of F/i-i-Z2, the SiH" (SuH") branch of simple (super) hexagons is on 
the group orbit of the SiH"*" (SuH"'") branch since k(z) = — z. Note that the only states that are 
the same (after rescaling x) for every value of a and (3 are the rolls, the rhombs Rh/io, the simple 
hexagons, and the simple triangles. Rolls are the only state that is common to both the square and 
the hexagonal lattices. 

Only the super- and the anti-states of Tables 7 and 8 are characterized by translation free 
isotropy subgroups. Hence these are the only axial planforms with (smallest) periodicity determined 
by ii and £2. All of the other axial planforms are periodic on a finer square, hexagonal, or rhombic 
lattice. In particular, the wavelength of their periodicity is 1/kc, where kc = \/ o? \ fP' for square 
lattice states and = \/ o? -\- f]"^ — a(3 for hexagonal lattice states. Note that while the periodicity 
of super- and anti-states is given by = \l2\ ^ l/^c, the lengthscale l/fcc is also evident in the 
patterns. This lengthscale shows up as small scale structure in the patterns; compare, for example, 
simple hexagons with super hexagons in Figure 6. 

The axial planforms, rhombs Rh/io, simple squares, and simple triangles/hexagons, are listed 
in Tables 3, 4 and 6 as possessing Z„ symmetry, where n = 2,4,6, respectively. However, on 
the appropriate (finer) square or hexagonal lattice that supports these planforms, they are D„- 
symmetric. Once we include the hidden refiection symmetry t^i we recover the full D„ symmetry 
that is manifest in Figures 5-7. 

Finally, we note that there is a countable set of rhombs that are periodic on square or hexagonal 
lattices. In Table 9 we characterize the rhombs on the square and hexagonal lattices in two different 
ways. We give the angle between the wave vectors associated with the critical modes, e.g. the angle 
between Ki and K3 for ^\^sl,a,|3■ We also give the aspect ratio of the rectangles evident in the 
rhomb patterns in Figures 5-7. 
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Table 8: Axia 


Planforms for 12- dimensional representations of Th anc 


Th + Z2. 


Name 


Nomenclature 


representative point z G Fix(S) C 


r 


Rolls 


R 


z = (1,0,0,0,0,0) 


Th and Th + Z2 


Simple Hexagons 


SiH± 


z = ±(l, 1,1, 0,0,0) 


Th and Th + Z2 


Simple Triangles 


SiT 


z = (i, i, i, 0, 0, 0) 


Th + Z2 only 


Rhombs 




z = (l, 1,0,0,0,0) 


T/, + Z2 only 


Rhombs 




z = (l, 0,0,1, 0,0) 


T/, and Th + Z2 


Rhombs 




z = (l, 0,0,0,1,0) 


T/, and Th + Z2 


Rhombs 




z = (l, 0,0,0,0,1) 


T/, and Th + Z2 


Super Hexagons (''^ 




z = ±(1,1, 1,1, 1,1) 


T/, and Th + Z2 


Anti-hexagons 


AH„^^ 


z = (l, 1,1, -1,-1,-1) 


T/, Z2 only 


Super Triangles 


SuT„^^ 


z = (i, i, i, i, i, i) 


Th + Z2 only 


Anti-triangles 


AT„,^ 


z = (i, i, i, —i, —i, —i) 


Th + Z2 only 



Simple hexagons are simply "hexagons" in [13]. 
Simple triangles are called "right triangles" in [13]. 
Golubitsky et al. call this rhombs state the "patchwork quilt" [13]. 
(''^ Super hexagons are called "hexagons" in Dionne and Golubitsky [7]. 



Table 9: Characterization of the rhombs. 



Lattice 


Rhombs 


aspect ratio 


angle 


Square 




a-13 
a+P 


(^) 


Square 


Rhs2,a,/3 


P 
a. 




Hexagonal 


^^hl,a,l3 


2fj-a 
\/S a 


/ a' + 2al3-20' \ 
\2(a^-al3+l3^) 1 


Hexagonal 


Rh/i2,a,/3 


a+l3 




Hexagonal 


^^h?,,a,l3 


2a-P 
V3 


COS-1 i 2«'-2«/3-/3^ \ 


Hexagonal 




1 

^/3 


7r 
3 
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4 Stability Results: Square Lattice. 



In this section we compute the linear stability, at bifurcation, of the six axial planforms listed in 
Table 7. We do this within the center manifold framework of a general Fs-equivariant bifurcation 
problem z = g(z,A), g : X R ^ C*. An equilibrium solution branch za of the bifurcation 
problem is linearly (orbitally) stable if all eigenvalues of the Jacobian matrix Dg(zA, A), not forced 
by symmetry to be zero, have negative real part for A sufficiently close to zero. If any eigenvalue 
has positive real part then the planform is unstable. 

We begin by determining the restrictions that symmetry imposes on the 8x8 (real) Jacobian 
matrix Dg evaluated on each of the axial solution branches. This allows us to determine the 
eigenvalues of the Jacobian matrix in terms of certain entries. We then determine the Taylor 
expansion, about z = 0, of a smooth r^-equivariant vector field to high enough order to determine 
the sign of the real part of each eigenvalue. 



4.1 Commuting linear maps. 

Let gj = g'j -\- ig'^^ j = 1, ...,4, where the r and i superscripts specify the real and imaginary parts, 
and let Zj = Xj + iyj, where Xj and yj are the real and imaginary parts of Zj, respectively. Thus 
Xj = g'j and i/j = g'j, j = 1, ...,4, are the eight components of the equivariant vector field over the 
reals. Table 10 gives the eigenvalues of the Jacobian matrix evaluated on each of the axial solution 
branches. We use two different approaches to determining the eigenvalues. 

The first approach exploits the observation that the Jacobian matrix evaluated on a solution 
branch z\ commutes with each element a G ^zx- For example, the Jacobian matrix Dg evaluated on 
the rolls solution branch commutes with the linear transformations associated with the generators 
of S = Z2-i-S^ in Table 3, namely 

z^z (4.1) 

and 

z ^ {zi,e^^'^"'+^'>Z2,e^^'^"'-^'>zs,e*^'"^'z4), s G R. (4.2) 

We choose an ordering of the coordinates of R® to be ixi,X2,X3,X4,yi,y2,y3,y4). It follows from 
the observation that Dg, evaluated on the rolls solution branch, must commute with the above 
transformations that Dg is diagonal and three of the eigenvalues have multiplicity two. Moreover, 
the group orbit of rolls is one- dimensional so there is a zero eigenvalue associated with translation 
along the group orbit. The null direction is determined by computing the tangent vector to the 
group orbit, e.g., 

d 

0(a;i, 0, 0, 0) = (-2Triaxi,0, 0, 0), G R , (4.3) 



801 



e=o 



where the action of on is given by (2.25). It follows that (0,0,0,0,1,0,0,0)''' is a null 

eigenvector of Dg and that the eigenvalue is zero. 

The second approach to computing the eigenvalues relies on forming the isotypic decomposition 
of for the isotropy subgroup of a solution z\. This decomposition determines coordinates 
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that block-diagonalize Dg [12]. The isotypic decomposition proceeds by first decomposing 
into S-irreducible subspaces Vj so that = Vq ® ® • • • ® V^. (Recall that a representation is 
S-irreducible if the only S-invariant subspace of Vj, other than {0}, is Vj itself.) The isotypic 
components Wj are then formed by combining the irreducible subspaces that are T,- isomorphic. 
(Two S-irreducible subspaces are S-isomorphic if there exists a linear isomorphic mapping between 
them which commutes with the action of S.) The isotypic decomposition is = l^offil^i®' ■ '®Wk, 
k < i, where the Wj are uniquely determined. 

As an example of the second approach, which block-diagonalizes Dg, we determine the isotypic 
components of for S = D4[i?^/2 7 ''"a^i] that applies to the super squares state. A decomposition 
of into irreducible subspaces is 

= R{(1,1,1,1)}®R{(1,1,-1,-1)}®R{(1,-1,1,-1)}®R{(1, -1,-1,1)} 

®R{(i,i,i,i),(-i,i,-i,i)}® R{(i,0,0,i),(0,i,-i,0)} . (4.4) 

The one- dimensional subspaces are non-isomorphic representations of D4, whereas the two-dimensional 
subspaces are D4-isomorphic and hence are in the same isotypic component (identify (i, i, i, i) with 
(i,0,0,i) and ( — with (0,i,— i,0)). From the one- dimensional isotypic components we 

can immediately determine four of the eigenvalues of Dg evaluated on the super squares solution 
branch; these are 

^ + ^ + ^ + ^ ^ffi I ^ffi ^ffi ^ffi 
dxi dx2 dxs 8x4 ' dxi 8x2 dx^ 8x4 ' 
dgl _ dgl ^ dgl _ dgl dgl _ dgl _ dgl ^ dgl 
dxi dx2 8x3 8x4 ' dxi 8x2 8x3 8x4 

Symmetry places additional restrictions on the matrix obtained by restricting Dg to the four- 
dimensional isotypic component W4. Specifically, Dg|^ commutes with the action of D4 on W4 = 
R"*, where W4 is the direct sum of two isomorphic D4-absolutely irreducible subspaces. (Recall 
that a representation of a group T acts absolutely irreducibly on a space V if the only linear maps 
on V commuting with T are multiples of the identity.) It follows that A = Dg|^ has the form 

where I2 is the 2x2 identity matrix and a, b,c,d £ R. Each eigenvalue in this matrix has multiplicity 
two. Moreover, two of the eigenvalues must be zero because the group orbit of super squares is 
two-dimensional. Thus the eigenvalues are determined by simply computing Tr(A), which, in terms 
of real coordinates, is 

Tr(A) = M + M + M + |^. (4.7) 
8yi 8y2 8y3 8y4 

This can be further simplified by noting that since Dg commutes with the transformations i?^/2 

dgl ^ dgl ^ dgl ^ dg\ 
dyi dy2 dy3 dy4 



and 
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on the super squares solution branch. Thus the repeated eigenvalue, ^Tr(A), is simply 2|^. 

The isotypic decomposition of for S = D4 is the same as the D4-isotypic decomposition 
presented above. The details of the computations of the eigenvalues for the remaining axial plan- 
forms are omitted. Note that symmetry considerations alone determine that the eigenvalues of Dg 
are real for all of the axial planforms. 

4.2 Normal form for D4-i-T^. 

In this section we determine the Taylor expansion of the equivariant bifurcation problem (2.16) to 
sufficient order to determine the signs of the eigenvalues of Dg given in Table 10. 

The equivariance condition (2.17) for F = is satisfied if (see, for example. Appendix A. 3 in 

[6]) 

Zl = 9l(zi^ Z3^ Z4) 
Z2 = 9i(z2t Zi^ Z4^ Z3) 

Z3 = gi(z3,^'i,Zi,J2) (4.9) 
Z4 = 9l(z4^ Z3, Z2, Zl) , 

where 



5i(z) = 5i(z), (4.10) 

and 

QUgiiz)) = zm{z) , for all G . (4.11) 



Equivariance with respect to D4 C F^ is guaranteed by conditions (4.9) and (4.10), while equivari- 
ance with respect to C F^ is equivalent to condition (4.11), i.e., to 'zigi(z) being an invariant 
function of the T^-action. The problem of determining the general form of the Fg-equivariant vector 
field then reduces to one of finding the most general function gi(z) that satisfies (4.10) and (4.11). 

We assume that local to the bifurcation point z = 0, A = the equivariant normal form can 
be expanded in a Taylor series about z = 0. We proceed by determining the general form of a 
T^-invariant function h = zigi(z). Let 

Hz) = ^ ak z'l^z'l-zl^4^ztztz',^4^ , (4.12) 

fc,>0 

where k2 > 0. The condition (4.10) determines that the coefficients U]^ are real. Given the action 
(2.25) of on C*, the T^-invariance of h determines that ay^ = unless 

ih - k2) ia0i + (302) + ih - h) (-(304 + a02) 
+ ih-ke)ip0i + a02) + ikr-ks)i-a0i+(302) = (4.13) 

for all (01,02) G T^. Clearly, ki = k2, k^ = k^, k^ = kg, and ki = kg is a solution to equation 
(4.13), which yields the T^-invariants \zj\'^, j = 1, 2, 3, 4. In the following we factor out all powers of 
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Table 10: Eigenvalues and their multiplicities for axial planforms associated with 8-dim. represen- 
tations of r,. 



Axial Planform 



Eigenvalues 



Rolls 

5(1,0,0,0), X eK 



M 

dxi ' dx2 



(mult. 2), 1$ (mult. 2), |$ (mult. 2), 



Simple Squares 
: a;(l, 1,0,0), X eK 

S = Z4-I-S1 2 



dxi ' dx2 ' dx2 ' 9a:3 



(mult. 4), (mult. 2) 



Rhombs (Rhsi,^,/?) 
= a;(l, 0,1,0), a; G R 
S = D^+Si,3 



dx\ dxs ' dxi dxs ' 9a:2 



(mult. 4), (mult. 2) 



Rhombs (Rhs2,a,/3) 
z = a;(l, 0,0,1), X eK 
S = Df +Si,4 



dxi 



9^:4 ' dxi 



M 

8X2 



(mult. 4), (mult. 2) 



Super Squares 
z = x(l, 1, 1, 1), a; G R 
S = D4 



9a[ 


, 9ar 


9a[ 


, 9ar 


9a[ 


, 9ar 


9a[ 


9a[ 








' 8xi ' 




' 9a;2 




9^:4 




_ M _ 


L M 


_ M 




_ M _ 


M _ 


L M 




8X2 


9a;3 


9^:4 ' 




8x2 


8x3 


9^:4 



2^ (mult. 2), (mult. 2) 



Anti- Squares 
:i, 1,-1,^-1), a; G R 
S = D4 



Same as Super Squares 
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I^TjP from the translation invariant monomials P]^(z) = z^^ z^^ Z2'^ z^^ z^^ ^ and only consider 
monomials of the form z^Z2Z^zl^ where we adopt the convention that zj^ = if m < 0. In this 
way we reduce the problem of finding all T^-invariant monomials to one of finding m,ra,p, g G Z 
such that 

m(a^i + /3^2) + n(-/3^i + 0^2) + + "^2) + + ^62) = 0. (4.14) 

Since 9i and 62 are independent this requires 

(m — q)a — (n — p)(3 = , (n + p)a + (m + g)/3 = . (4-15) 

Furthermore, since (a, (3) = 1, this implies 

m — q = j (3 , n — p = ja , (4.16) 

n + p = k(3 , m + q = —ka , 

where j, k £ Z. Solving for m, n,p and q gives 

m= -^(A;a-i/3) , n = ^(ja + A;/3) , (4.17) 

P = -^(i" - 5 ? = -^(^a + i/3), 

Since (q!,/3) = 1 and a and (3 are not both odd, the system of equations (4.17) has no (nontrivial) 
solution if any two of m, n,p or q are zero. It then follows that both k and j are even and equations 
(4.17) may be replaced by 

m = -{k'a - j'fi) , n = (j'a + k' fi) , (4.18) 

P =-(/«- A;'/3) , q = -{k'a+j'(3), 

where j' ,k' G Z. If j' = 0,A;' 7^ or j' ^ 0,A;' = then (4.18) yields the translation invariant 
monomials 

^1 ^2 % '^4 ' ^1 ^2 '^3 '^4 ' (4.19) 

and their complex conjugates. 

The monomials (4.19) are order 2{a-\- (i). We now show that these are the lowest order (nontriv- 
ial) translation invariant monomials of the form z'^Z2Z^z\. To do this we consider the remaining 
cases for which j'k' 7^ in 4.18: 

1. j' > 0,A;' > 

2. j' > 0,A;' < 

3. j' < 0,A;' > 

4. j' < 0,A;' < 0. 
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In case 1, since no solution to (4.18) exists if two of m, n,p and q are zero, the order of the invariant 
is 

\m\ + \n\ + IpI + |g| > \n\ + |g| 

> (f + k')ia + (3) 

> 2{a + l3) since j',k'>0, 

Similarly, in case 2 the order of the invariant 

\m\ + \n\ + IpI + |g| > \m\ + |p| 

> (j'+\k'\)(a + f^) 

> 2(a + /3). 

The argument that \m\ + \n\ + |p| + |g| > 2(a + (3) is similar for cases 3 and 4. Hence no invariants 
of order less than or equal to 2(a + (3) occur for j'k' 0. 

We use (4.11) to compute the leading order terms in the Taylor expansion of gi(z) from the 
T^-invariants, \zj\'^, (4.19), and their complex conjugates. In particular, we find 

ii = ^i/(kiP, \Z2\\ \Z3\\ \Z4\^) + b^4-^z^z^z^^ + b^zr^z^^z^zl + 0(2(a + /3)), (4.20) 

where it follows from (4.10) that / is a real-valued function of its arguments and that 61,62 G 
Condition (4.9) determines the remaining components of g from gi. 



4.3 Stability. 



We use the leading order terms in the Taylor expansion of the normal form (4.20) and the expressions 
for the eigenvalues given in table 10 to compute the signs of the eigenvalues of Dg at bifurcation. 
The results for the axial planforms are summarized in Table 11. The eigenvalues for the rolls, 
simple squares and rhombs are determined using a cubic truncation of the normal form. 



ii = A2;i + 2;i(ai|2;ip + a2|2;2p + asksP + a4k4p) + C'(lzl^) 

Z2 = Xz2 + Z2iai\z2\'^ + a2\zi\'^ + aslzil"^ + ttilzsl"^) + 0(\zf] 

Z3 = Xz3 + zsiailzsl"^ + a2\z4\'^ + aslzil"^ + a4\z2\'^) + 0(\z\^] 

Z4 = \z4 + Z4{ai\z4\'^ + a2\z3\'^ + a3\z2\'^ + a^lzil"^) + 0{\z\^] 



(4.21) 



Here we assumed, without loss of generality, that time has been scaled so that the linear term in 
the normal form is simply Az. We note that the computation of the last (repeated) eigenvalue for 
the super squares and anti-squares solution branches is simplified by the observation that 



dgi _ dgi 

dzj dzj 



(4.22) 



The sign of this eigenvalue is determined by keeping all terms through 0(2(a + (3) — l) in the Taylor 
expansion of the normal form (4.20). 
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Table 11: Stability results for the square lattice bifurcation problem, from Table 10, and equations 
4.20, 4.21. 



Axial Planform 


Signs of Nonzero Eigenvalues 


Branching Equation 


Rolls 
z = (s, 0,0,0) 


sgn(ai), sgn(a2 — ai), 
sgn(a3 — ai), sgn(a4 — ai) 


= Aa; + ais"^ + • • • 


Simple Squares 
z = (x, X, 0, 0) 


sgn(ai + 02), sgn(ai — 02), 
sgn(a3 + 04 — ai — 02) 


= Aa; + (ai + 02)2;"^ + • • • 


Rhombs (Rhsi,^,/?) 
z = (x, 0, X, 0) 


sgn(ai + 03), sgn(ai — 03), 
sgn(a2 + 04 — ai — 03) 


= Aa; + (ai + 03)2;"^ + • • • 


Rhombs (Rhs2,a,/3) 
z = (x, 0, 0, x) 


sgn(ai + 04), sgn(ai — 04), 
sgn(a2 + — ai — 04) 


= Aa; + (ai + 04)2;"^ + • • • 


Super Squares 


sgn(ai + a2 + + 04), 
sgn(ai + a2 — — 04), 
sgn(ai + 03 — 02 — a4), 
sgn(ai + 04 — 02 — a3), 
sgn( — bi(3 — b2a) 


= Aa; + (ai + a2 + f^s + a4)a;"^ + • • • 
+(5i + 52)a;2("+'^)-^ + ••• 


Anti-squares 

2 — ^ tZ/ ^ tZ/ ^ tZ/ ^ tZ/ ^ 


sgn(ai + a2 + + 04), 
sgn(ai + a2 — — 04), 
sgn(ai + 03 — 02 — a4), 
sgn(ai + 04 — 02 — a3), 
sgn(bi(3 + 52Ci) 


= Aa; + (ai + a2 + f^s + a4)a;"^ + • • • 
-(hi + b2)x^("+^'}-^ + ■■■ 
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We assume that the following nondegeneracy conditions are satisfied: 



«! 7^ 0, ±a25 ±«35 ±«4 5 

(ai + 02) 7^ ±(a3 + ^4) , 



5i/3 + 7^ , 



(4.23) 



In this case we can draw a number of conclusions from Table 11. 

1. Any one of the axial solution branches can bifurcate supercritically to produce a stable solu- 
tion. 

2. If the super squares and anti-squares are neutrally stable at cubic order, then one and only 
one of these two states bifurcates stably. 

3. If all of the axial planforms bifurcate supercritically, then at least one of them must be stable. 

4. If any axial solution branch bifurcates sub critically, then rolls, super squares and anti-squares 
are all unstable. 

5. If rolls, super squares, or anti-squares bifurcate sub critically, then all axial planforms are 
unstable at bifurcation. 

6. If simple squares is the only axial solution branch to bifurcate sub critically, then it is still 
possible that one, but not both, of the rhombs solutions is stable. Similarly, if one of the 
rhombs solutions bifurcates sub critically, then it is possible that simple squares or the other 
rhombs solution branch is stable, though they cannot both be stable in this case. 

7. The only solution branches that can co-exist stably are simple squares SiS and the rhombs 
^^si,a,fj, Rhs2,a,/3- Any combination of two of these states can bifurcate stably, but not all 



5 Stability Results: Hexagonal Lattice. 

In this section we compute the linear stability, at bifurcation, of the axial planforms that are 
associated with the twelve- dimensional representations of Th and Th + Z2 (see Table 8). As with 
the square lattice case, we do this within the framework of a general F-equivariant bifurcation 
problem z = g(z. A), where g : X R ^ C^. 

In the case of T = Th + Z2 there are only odd terms in the Taylor expansion of g due to the Z2 
symmetry. However, if the Z2 symmetry is absent, then even terms are admissible. In particular, 
we find that the coefficients of most, but not all, quadratic terms in the Taylor expansion of gi are 
zero; the exception is e = 2 8^2 8^ ' following vector is F/j-equivariant 



three. 
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The presence of such a quadratic term ensures that generically all of the axial planforms bifurcate 
unstably [15]. In order to obtain stable axial solution branches we focus on the degenerate bifur- 
cation problem defined by e = 0. We then discuss briefiy the unfolding of this bifurcation problem 
(i.e., the case < |e| ^ 1), before analyzing the generic Th + Z2-equivariant bifurcation problem. 



5.1 Commuting linear maps. 

We begin by determining the restrictions that symmetry places on the eigenvalues of the 12 X 12 real 
Jacobian matrix Dg when it is evaluated on an axial solution branch. The results are summarized 
in Table 12. We denote the real and imaginary parts of gj by g'j and g'j, respectively, and the real 
and imaginary parts of Zj by Xj and yj, respectively, j = 1, ...,6. 

We exploit the observation that Dg, evaluated at z, commutes with the generators of the 
isotropy subgroup Sz to determine the general form of the Jacobian matrix. For example, in the 
case of rhombs, Khhi^a,fj, we find that the Jacobian matrix has the form 



Dg 



A 
B 



(5.2) 



where the 6x6 real matrices A and B have the form 



B 



( (111 





ai4 







\ 







(122 


a23 















(132 


^33 












(ll4 





an 










•) 











^33 


032 




^ 








O23 


022/ 












5l4 







^ 





(122 


-0.23 
















-0.32 


O33 













5l4 








hii 
























^33 


-032 


^ 













-023 


022 I 



(5.3) 



Here we have chosen an ordering of the coordinates of R^^ given by (xi, a:2, ^3, ^4, ^5, ^e? Vx-) 
2/2, 2/3, 2/4, 2/5, ye), so, for example, an = and a23 = g^- Moreover, by determining the 
null directions associated with the two-dimensional group orbit of Khhi^a/J solutions we find that 
bii = bi4 = 0. The computations for the other axial solution branches are similar. Note that rolls 
and simple hexagons lie in the six- dimensional fixed-point subspace Fix(Si^2,3) on which the hidden 
refiection t^^ acts (equation 2.28); for these solutions we take the hidden symmetry into account 
in determining the eigenvalues of Dg. 

In the case of super hexagons, we find the computation of the eigenvalues of Dg is simplified 
by forming the De-isotypic decomposition of C^. It is 



R(l,l,l,l,l,l)® R(l, 1,1, -1,-1,-1) ®R(i,i,i,i, i,i)®R(i,i, i,-i,-i,-i) 
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® R{(1, -1,0,0, -1,1), (-1,0, 1,-1, 1,0), (0,1, -1,-1, 1,0), (1,-1, 0,1,0,-1)} 

® R{(-i,0,i,i,0,-i),(0,-i,i,0,i,-i),(0,i,-i,-i,i,0),(-i,i,0,-i,0,i)} . (5.4) 

The 12 X 12 Jacobian matrix Dg is block diagonal with respect to this basis [12]. The four eigen- 
values associated with the four one- dimensional isotypic components are thereby determined. In 
addition, each four- dimensional isotypic component is a direct sum of two isomorphic De-absolutely 
irreducible subspaces of dimension two. Thus the Jacobian matrix, restricted to either of these sub- 
spaces, must have the form 

ah bl2 



where I2 is the 2x2 identity matrix, and a,b,c,d G R. The eigenvalues, 121,122, each have multi- 
plicity two and satisfy = ad — he, jii -\- 112 = a + d. Moreover, one of the repeated eigenvalues 
associated with the last isotypic component in (5.4) must be zero due to the symmetry; thus 
the nonzero eigenvalue (of multiplicity two) is determined by computing the trace of Dg restricted 
to this isotypic component. 

5.2 Normal Form for Dg+T^. 

In this section we determine the Taylor expansion of the equivariant bifurcation problem (2.16) to 
sufficient order to determine the signs of the eigenvalues of Dg given in Table 12. We focus, in 
particular, on the degenerate bifurcation problem e = 0, where e is the coefficient of the quadratic 
term (5.1). 

Our approach to determining the leading terms in the F/j-equivariant normal form is the same 
as that employed in Section 4.3 for the Fg-equivariant normal form. The general F/^-equivariant 
vector field that satisfies the equivariance condition (2.17) is 

Zl = 9l(zi, Z2, Z3, Z4, Z5, Zq) 

Z2 = 9i(z2t Z3, Zl, Z5, Zq, Z4) 

Z3 = 9i(z3t Zl, Z2, Zq, Z4, Z5) 

Za = gi(z4,Z6,Z5,Zi,Z3,Z2) (5.6) 

Z5 = 9i(z5t Z4, Zq, Z2, Zl, Z3) 

Z& = gi{zQ, Z^, Z4, Z3, Z2, Zl) , 

where 
and 



gi{7,) = gi{z) , (5.7) 

<d(zigi{z)) = zigiiz) , for all G . (5.8) 

Equivariance with respect to De C F/^ is guaranteed by conditions (5.6) and (5.7), while equivariance 
with respect to C F/^, with action given by (2.32), is equivalent to condition (5.8). Finally, the 
hidden refiection (3.1) puts an additional restriction on the function gi(z), namely 

gi(zi,Z2, Z3, 0, 0, 0) = gi(zi,Z3, Z2, 0, 0, 0) (5.9) 
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Table 12: Eigenvalues for axial planforms associated with 12-dim. representations of F/j. 



Axial Planform 



Eigenvalues 



Rolls 

z = a;(l, 0,0, 0,0,0) 



'^'^ (mult. 2), 1$ (mult. 2), |$ (mult. 2), 



M + M (a) (mult. 2), 1^ 



M (a) 



dxe 

(mult. 2), 



Simple Hexagons 
z = a;(l, 1,1, 0,0,0) 

S = Ze+Si 2,3 



M + 2M(»), M - M (») (mult. 2), 



If (mult. 6), 



(mult. 2) 



Rhombs (Rh/,i,„,^) 
z = a;(l, 0,0, 1,0,0) 
S = D-+Si,4 



dx^ 



(mult. 2) 



//i,//2 (mult. 4); //I + = U + If , = Ulf - glf 



Rhombs (Rh/,2,a,/3) 
z = a;(l, 0,0, 0,1,0) 
S = D^+Si,5 



Iii,li2 (mult. 4); //I + //2 



M 

dx2 



M 

9a;3 



(mult. 2 



9g; 9fl; 

9a;2 9a;3 



9g; 9fl; 
9a;3 9a;2 



Rhombs (Rh/,3,„,^) 
z = a;(l, 0,0, 0,0,1) 
S = Df +Si,6 



9a;6 



//i,//2 (mult. 4); //I + 112 



M _ 
_ M 

9a;2 



9a;6 



(mult. 2) 



9a5 

9^ ' /^l/^2 



8^2 9^3 
9a;2 9a;3 



8^2 9^3 
9a;3 9a;2 





_ M 


_ M _| 


_ M J 


L M _ 


L M 


M _ 


L M _ 


1 M 


_ M 


_ M 


_ M 


dx\ 


9a;2 


9a;3 






I" 9a;6 ' 




9a;2 


9a;3 




9a;5 


9a;6 


dg\ 


J_ M . 


J_ M . 


j_ M _ 


, M _ 


1 M 


M _ 


lM , 


_ M 


_ M 


_ M 


_ M 


dyi 


9|/2 


9|/3 


9|/4 


9|/5 


^ 9|/6 ' 


dyi 


^ 9|/2 ^ 


dys 


dyi 


dys 


dye 



Super Hexagons 
z = a;(l, 1,1, 1,1,1) 
S = De 



' 9|/i 



M _ M (mult. 2), 

9|/2 9|/3 V 



(mult. 2), 1^1, (mult. 2) 



9a:i 



9a;2 , 



9a:i 



9a;3 



8X2 



9a;3 , 



dx^ 



9a;5 



9^:4 dxQ J \.dxs dxQ 



-) } , Ml + 

Here the effect on Dg of the hidden symmetry t^^ (3-1) is included. 



9:ri 9:^2 9:^3 
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The problem of determining the general form of the F/i-equivariant vector field then reduces to one 
of finding the function gi(z) that satisfies equations (5.7)-(5.9). 

As in the square lattice case, we assume that local to the bifurcation point we can Taylor expand 
the function gi(z). We proceed by determining the T^-invariant monomials, which represent all 
terms present in the Taylor expansion of the T^-invariant function zigi(z). In the following, we 
assume that overall factors of \zj\'^, j = 1,...,6, which are manifestly T^-invariant, have been 
removed from the monomials that we consider. We follow the convention that z'j = Zj if n < 0, 
and focus on T^-invariants of the form 



"2 ^3 ^4 ^5 ^6 5 

where m, n,p,q,s £ Z satisfy 



zV^zliz^^zUlzl, (5.10) 



(m — n + q — s)a + (n — p — r + s)[i = 0, 
— {n—p—q-\-r)a-\-{m — p—q-\-s)[i = 0. (5-11) 

Since a and [i are relatively prime, we have 

m — n + q — s=j(3, n — p— r + s = —ja , (5.12) 

n — p — q + r = k(3, m — p — q + s = ka , 

where j,k G Z. 

There are no nontrivial solutions of (5.12) with more than three of m,n,p),q,r, s zero. In the 
case that j = A; = in (5.12), then m = n = p and q = r = s, which yield the invariants 

2:12:2-2^3 and Z4Z5Ze , (5.13) 

and their complex conjugates. 

We claim that the lowest order (a, /3)-dependent invariants + \k\ 0) are 

J3 -a-l3 -a-l3 13 -,l3 -a - 13 -a - 13 13 -,l3 -a - 13 -a - 13 13 (K 1A\ 

^1^2 ^4 ^5 5 ^2^3 ^5 ^6 ' ^3^1 ^6 ^4 ' W"-'-^/' 

and their complex conjugates, which are order 2a. Note that the set of the three invariants (5.14) 
together with their complex conjugates is invariant under the action of De. We justify the above 
assertion by the following steps: 

1. Show that we can assume that one of m, n,p is zero and that one of q, r, s is zero. We then 
focus on the specific case with p = since the invariants with m = and ra = can be 
transformed to p = by the action of i?^/3 G De on the invariant. 

2. Consider the cases where p = and exactly two of m, ra, q, r, s are zero. 

3. Consider the cases where p = and only one of q, r, s is zero. 
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Step 1. Consider m,n,p > 0. If m,n,p are all nonzero, then we can construct a lower degree 
T^-invariant monomial from (5.10) by factoring out the invariant ziZ2Z^. Thus the lowest degree 
monomials cannot have m,n,p > 0. Similarly, if m,n,p < 0, then we can lower the degree by 
factoring out the invariant ziZ2Z'i. 

Suppose now that m, ra,p do not all have the same sign; for example, consider the case m > n > 
> p. Then the order of the monomial is m + ra + |p| + |g| + \r \ + However, if z^zl^^zj^^zlz^Zg is 
invariant, then which is of lower order m + |p| + |g| + \r\ + |s| unless ra = 0, 

in which case it has the same order. Since we only aim to find the lowest order (a, /3)-dependent 
monomials, we can assume ra = in (5.10). The argument is the same for the other orderings of 
m,ra,p, and 0; in each case we can find a lower degree invariant monomial unless one of m,ra,p is 
zero. 

Hence, the lowest order (a, /3)-dependent monomial has at least one of m, ra, p equal to zero, and 
by a similar argument we can assume that one of q,r,s is zero. In the following steps, we assume 
that p = 0. The invariants with m = or ra = and p are obtained from the invariants with 
p = by applying i?^/3 G De to the monomials. 

Step 2. Let p = and exactly two of m,ra,g,r, s be zero. From step 1, we can assume that 
at least one of q, r, s is zero. There are nine combinations to consider; in each case we obtain an 
invariant of degree greater than 2a. As an example, consider p = r = s = 0. Then equations (5.12) 
can be solved provided we choose j and k such that 

(k-j)a = (j + 2k)f:S. (5.15) 

Thus 

k-j = l/3, j + 2k = la, (5.16) 
where / G Z. Solving for j and k gives 

j = h(a-2/3), k = h(a + /3), (5.17) 

and, since (3, a + /3) = 1, / must be divisible by 3. Let / = 3/', /' G Z, then we find 

m = l'fi(2a - fi) , 

ra = l'a(2fi-a), (5.18) 

q = -/'(q,2 _ Q,^ + . 

The order of the invariant z^Z2zl is 3l'a(3, which is greater than the order of the invariants (5.14). 
The other eight combinations, with three non-zero exponents, also give invariants of order either 
3a(3, + a/3, or 2q!2 -|- a(3 — (3'^, all of which are greater than 2a. 

Step 3. If p = then it follows from (5.12) that 

m = j/3 - ^(j - k){a + (i) , 
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3 

We set (k — j) = 3/, / G Z, because (3, a + (3) = 1. Hence, 



kf:S-ja + ^(j-k)(a + f^), (5.19) 
j/3 - ka- ^{j - k){a + f3) , 
kf3 + lij-k)ia + (3) . 



m = + l{a + /3) , 
n = -j(a - /3) + 1(2/3 - a) , 
q-s = -j(a- fi)-l(2a- fi) , (5.20) 
r - s = j/3 + 1(2/3 - a) , 

where j, / G Z. 

We consider the three cases g = 0, r = 0, and s = separately. The invariant z^Z2~^z'^~^z^ 
in (5.14) is obtained in the case p = s = for / = 0, j = 1 in (5.20). The cases p = q = and 
p = r = lead to nontrivial invariant monomials of degree greater than 2a. 

For p = g = 0, the degree of the monomial is |m| + |ra| + |r| + |s|. It follows from the restriction 
a > /3 > a/2 in Table 2 that 

\m\ + \n\ + \r\ + \s\ > \m + n + r\ = \l\(a + 4/3) > 2\l\a . (5.21) 

Hence, the degree of the monomial is greater than 2a unless / = 0. However, if / = 0, then 

\m\ + \n\ + \r\ + \s\ = \j\(3a - /3) > 2\j\a . (5.22) 

This proves that the monomials associated with p = q = 0, ninrs 0, are all of degree greater 
than 2a. The argument in the case p = r = is similar. 

From the invariant function zigi, we can compute the general form of the equivariant vector 
field through 0(2a — 1). Specifically, 

h = Zi /i(Mi,M2,M3,M4,M5,M6,gi,gi,g4,g4) + ^2^3 /2('Wl,M2,M3,M4,M5,M6,gi,gi,g4,g4) 

+ ei zr^'-^lzlzr^' + 62 zt^zr^'zr^'4 + 0(2a) , (5.23) 

where 

""j = 5 ft = ^1^22^3 , ^4 = Z4^Z5Zq , (5.24) 

and 64,62 G R are constants. It follows from (5.7) and (5.9), respectively that 



fj(ui,U2,U3,U4,U5,ue,qi,qi,q4,q4) = /j(mi, M2, Ms, 'W4, ""s, Me, ft, ft, ^4, ^4) , (5.25) 
/j(Mi,M2,M3,0,0,0,gi,gi,0,0) = /j(mi, M3, M2, 0, 0, 0, gi, gi, 0, 0) , j = 1,2 . 

We use the De symmetry to determine the other components of the normal form from (5.23), as 
indicated by (5.6). 
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Table 13: Stability results for the hexagonal lattice bifurcation problem in the degenerate case 
e = 0, from Table 12, and equations 5.6, 5.23 and 5.26. The branching equations are given in 
Table 14. 



Axial Planform 


Signs of Nonzero Eigenvalues 


Rolls 

z = 0,0, 0,0,0) 


sgn(ai), sgn(a4 — ai), sgn{a^ — ai), 
sgn(ae — ai), sgn(a2 — ai) 


Simple Hexagons 
z = x{l, 1, 1,0,0,0) 


sgn(ai + 2a2), sgn(ai — 02), 
sgn(a4 + + uq — ai — 2a2), sgn[x(ci — bi — 2b2)\ 


Rhombs (Khhi^a/j) 
z = 0,0, 1,0,0) 


sgn(ai + 04), sgn(ai — 04), 
sgn(a2 + — ai — 04), sgn(a2 + uq — ai — 04) 


Rhombs iKhh2,a,fj) 
z = 0,0, 0,1,0) 


sgn(ai + as), sgn(ai — as), 
sgn(a2 + a4 — ai — as), sgn(a2 + ae — ai — as) 


Rhombs iKhh3,a,fj) 
z = 0,0, 0,0,1) 


sgn(ai + ae), sgn(ai — ae), 
sgn(a2 + a4 — ai — ae), sgn(a2 + as — ai — ae) 


Super Hexagons 
z = a;(l, 1,1, 1,1,1) 


sgn(ai + 2a2 + a4 + as + ae), sgn(ai + 2a2 — a4 — as — ae) 

- 262 - 64 - ^5 - ^6 + ci + C2 - C3)], 

- 262 - 64 - ^5 - ^6 + ci - C2 + C3)], 
s£to[-(2q! - /3)ei - (a + /3)e2], 
sgn(pi + 112) = sgn(ai - a2), 

sgn{jiiii2) = s£'ra[2(ai — a2)^ — (a4 — as)^ — (a4 — a^)^ — (as — ae)^] 



5.3 Stability for degenerate bifurcation problem. 

Expanding (5.23) through quartic order, we obtain 

ii = \zi + 6^2-^3 + -2i(aikiP + a2k2p + a2k3p + a4k4p + asksP + aekeP) 

+ ^2^3(&lklP + ^2^2^ + ^2^3^ + ^4^4^ + ^5^5^ + ^6^6^) 
+ (Ci 2:1 2:2 2^3 + £22:4^:52:6 + C3^4^s^6) + C'(lzl^) . 



(5.26) 



(5.27) 



The remaining components of z = g(z) are determined from (5.26) using (5.6). Again we assume 
time has been scaled so that the linear term in g(z) is Az. Table 13 gives the signs of the eigenvalues 
in Table 12 for the degenerate case e = 0. The quartic truncation is sufficient to determine the sign 
of all eigenvalues except 

dg{ _ }_(dgj_ ^ dg\- 
dyi 2V%2 %3' 

for super hexagons. The sign of this eigenvalue is determined by retaining the leading order (q!,/3)- 
dependent terms in the normal form (5.23). 

We assume that e = and that the following nondegeneracy conditions are satisfied: 

ai 7^ 0, a2, ±a4, ±as, iae , 
(ai + 2a2) 7^ 0, ±(a4 + as + a^) , 
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2(ai - 
ci - bi 



ai — a2 




±(a4 - 05), ±(a4 - ae), ±(a5 
(a4 - a^f + (04 - ae)^ + (as - 
0, 64 + 65 + 56 ± (C2 - C3) , 



(5.28) 



a + /3 
2q! - /3 ■ 



In this case we can draw a number of conclusions from Table 13. 

1. While all axial solution branches bifurcate unstably when e 7^ 0, we find, in the degenerate 
case e = 0, that any one of the axial solution branches can bifurcate supercritically to produce 
a stable solution. 

2. There are two distinct branches of simple and super hexagons, denoted SiH^ and SuH^, 
respectively, associated with a; > and a; < 0. If simple hexagons are neutrally stable at 
cubic order, then one and only one of the two branches SiH^ is stable. If super hexagons are 
neutrally stable at cubic order, then one and only one of the two branches will be stable if 
(2q! — /3)ei + (a + /3)e2 > 0, while they are both unstable if (2q! — /3)ei + (a + /3)e2 < 0. 

3. If (2q! — /3)ei + (a + /3)e2) < then it is possible for all of the axial planforms to bifurcate 
supercritically, but none be stable. On the other hand, if (2q! — /3)ei + (a + /3)e2 > and all 
axial planforms bifurcate supercritically, then at least one of them must be stable. 

4. If any axial solution branch bifurcates sub critically, then rolls and super hexagons are unsta- 



5. If rolls or super hexagons bifurcate sub critically, then all axial planforms are unstable at 
bifurcation. 

6. If simple hexagons is the only axial solution branch to bifurcate sub critically, then it is still 
possible that one, but not more, of the rhombs solutions is stable. Similarly, if rhombs is the 
only axial solution branch to bifurcate sub critically, then it is possible for simple hexagons 
to be stable, or for one or more of the remaining rhombs solutions to be stable. However, if 
simple hexagons undone of the rhombs bifurcate sub critically, then all axial solution branches 
are unstable. 

7. If two of the rhombs solution branches bifurcate sub critically, then it is possible that the 
remaining rhombs solution or simple hexagons is stable, but not both. However, if all three 
rhombs solution branches are subcritical, then all axial planforms are unstable. 

8. The only solution branches that can co-exist stably are simple hexagons SiH and the rhombs 
^^hi,a,fj, Rh/i2,a,/3, ^^h3,a,fj- Any Combination of two of these states can bifurcate stably. It 
is also possible for all three types of rhombs to be stable simultaneously. However, if two or 
more of the rhombs are stable, then simple hexagons are unstable. 



ble. 
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5.4 Secondary Bifurcations. 

In this section we address briefly the unfolding of the degenerate bifurcation problem e = analyzed 
in the previous section. Speciflcally, we indicate how the stability of the axial solutions change along 
the solution branch in the case that |e| ^ 1. While a complete analysis of the unfolding is beyond 
the scope of the present paper, we do present an example in which part of a bifurcation diagram 
is computed. This example indicates the wealth of secondary transitions that occur close to A = 
when |e| ^ 1. When e 7^ 0, certain eigenvalues given in Table 13 are modifled to those given in 
Table 14. Note that, as discussed above, the presence of the quadratic term in the bifurcation 
problem ensures that at least one of the eigenvalues for each axial planform is positive for (A,z) 
sufficiently close to the origin. 

As a speciflc example, we consider the bifurcation problem 

zi = X zi + e Z2Z3 + zi (ail^^ip + a2k2p + a2k3p + a4k4p + as^sP + ae^eP) 

+ h Z2Z3 {\z2f + Izsf) + ei z^~'^~^z§z'^Zq~'^ + €2 z'^~'^ z!^''^ z'^~'^ z§ , (5.29) 

where 

ai = —1.5, a2 = —3.5, = 0.5, = 0.6, = 0.7, 62 = 0.6, ei = 1.0, 62 = 0.5. (5.30) 

This choice satisfles the non-degeneracy conditions (5.28). It follows from Table 13 that all three 
rhomb states are stable when e = 0. We show two bifurcation diagrams for < e ^ 1. The bifurca- 
tion diagrams indicate, schematically, the amplitude |z| as a function of the bifurcation parameter 
A for each axial planform. Solutions on the same group orbit are identifled and bifurcation points 
are indicated by solid circles. We follow the convention that solid lines indicate stable solutions and 
dotted lines indicate unstable solutions. Figure 8 is a well-known bifurcation diagram that applies 
to the six- dimensional representation of Th', here it is obtained by restricting our analysis to the 
six- dimensional subspace where z = (2:1,2:27-2^370,0,0). Figure 9 gives the bifurcation diagram that 
applies, for the same coefficient values (5.30), in the full twelve- dimensional space. 

In the six- dimensional subspace, where Z4 = Z5 = zq = 0, only two axial planforms exist, rolls 
and simple hexagons. In this subspace, and for the choice of coefficients (5.30), Figure 8 indicates 
that as A increases through 0, the trivial solution becomes unstable and there is a transition to stable 
simple hexagons. On further increase of A the hexagons become unstable and there is a transition to 
rolls. Both the transition to hexagons and that to rolls exhibit hysteresis. The bifurcation scenario 
of Figure 8 has been investigated in a wide variety of hydrodynamic systems [3, 10, 18, 23], in 
solidiflcation problems [28, 2, 21], and in chemical reaction- diffusion systems [9, 20]. 

Figure 9 indicates how the familiar bifurcation diagram in Figure 8 is modifled when we consider 
stability within the full twelve- dimensional space. In this case, rolls are always unstable to rhombs, 
and the range of stability of simple hexagons is greatly decreased. Indeed simple hexagons are 
stable only in a subcritical regime where super hexagons are also stable. In this case, on increasing 
A, there is flrst a jump at A = to stable super hexagons and then a transition to one of the three 
stable rhombs states. All of the transitions exhibit hysteresis. While the bifurcations to simple and 
super hexagons are transcritical, all other primary bifurcations are pitchforks. All of the secondary 
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Table 14: Stability results for the hexagonal lattice bifurcation problem in the case |e| ^ 1, from 
Table 12, and equations 5.6, 5.23 and 5.26. Also see Table 13; only the eigenvalues that depend on 
e are given here. 



Planform 


e-Dependent Eigenvalues 


Branching Equation 


Rolls 


ex + (a2 — ai)x^ + • • • 
— ex + (a2 — ai)x^ + • • • 


= Aa; + a4a;"^ 
+ C(a;5) 




ex + 2(ai + 2a2)a;^ + • • • 
— 2ex + 2(ai — 02)2;^ + • • • 
— ea; + (04 + + gq — ai — 2a2)x^ + • • • 
-3ea; + 3(ci - hi- 2h2)x^ + • • • 


= Aa; + ex"^ 
+ (a4 + 2a2)a;3 
+ C(a;4) 




Mi,M2; Ml + = (-2ai - 2a4 + 2a2 + as + ae)^;^ + • • • 
^4^2 = -e^x"^ 

-\-{ai + — a2 — as)('^i + — (12 — aejx'^ + • • • 


= Aa; + (a4 + a4)a;"^ 
+ C(a;5) 




ea; + 2(ai + 2a2 + (24 + as + ae)^^ + • • • 
ea; + 2(a4 + 2a2 — a4 — as — asjx^ + • • • 
-3[ea; + (64 + 2^2 + 64 + 5s + &6 - C4 - C2 + C3)x^] + ■■■ 
-3[ex + (bi + 21)2 + &4 + &5 + &6 - ci + C2 - C3)a;3] + • • • 
Mi,M2; Ml + = -4ea; + 4(a4 - a2)a;2 + • • • 
^4^2 = 4e^a;^ — 8(a4 — a2)ea;"^ + 4(a4 — a2)^x'^ 

-2[(a4 - as)^ + (a4 - a^f + (as - aQf]x'^ + • • • 


= Aa; + ex'^ 
+ (a4 + 2a2)a;3 
+ (a4 + as + a6)a;3 
+ C(a;4) 



(a) ^j^g results for Rh/i2,a,/3 (Rh/i3^„^^) are obtained from those for ^^hi,a,i3 by interchanging 



the 4 and 5 (4 and 6) subscripts. 

bifurcation points indicated in the diagram approach A = |z| = as e ^ 0. The paths of the 
secondary branches have not been computed. 

5.5 Stability results: F = F/j + Z2. 

In this section we consider the consequences of the additional Z2 symmetry, k(z) = — z, for the 
generic bifurcation problem on the hexagonal lattice, z = g(z,A), g : X R ^ C^. Specifically 
we consider the branching and stability assignments for the axial planforms listed in Table 8. 

The Z2 symmetry places some additional restrictions on the eigenvalues of rolls and rhombs 
listed in Table 12; specifically, it ensures that = = on these solution branches. The 
eigenvalues for simple and super hexagons, listed in Table 12, are unchanged. The eigenvalues of 
Dg for the remaining axial planforms are listed in Table 15. 

We note that the De-isotypic decomposition of is the same for the super hexagons, anti- 
hexagons, super-triangles and anti-triangles planforms; it is given by (5.4). Indeed, the only dif- 
ference between the eigenvalue structure for the triangle states and the hexagon states is that the 
null vectors lie in different isotypic components in the two cases. For example, the null vectors, 
associated with translations of the super hexagons and anti-hexagons, lie in the four- dimensional 
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isotypic component 



R{(-i,0,i,i,0,-i),(0,-i,i,0,i,-i),(0,i,-i,-i,i!,0),(-i!,i!, 0,-1!, 0,1!)} , 



(5.31) 



while the null vectors for the super and anti-triangles lie in 



R{(1, -1,0,0, -1,1), (-1,0, 1,-1, 1,0), (0,1, -1,-1, 1,0), (1,-1, 0,1, 0,-1)} . (5.32) 



The additional Z2 symmetry forces the coefficients of all even order terms in the Taylor expan- 
sion of the F/j-equivariant normal form (5.23) to be zero. Hence there are no quadratic terms; the 
differences between the degenerate bifurcation problem with F/j-symmetry (e = 0) and the generic 
bifurcation problem with Th + Z2-symmetry arise at (9(|z|'*). Thus the eigenvalues for the rolls 
and the rhombs in Table 13, which are determined by a cubic truncation of the normal form, are 
unchanged by the extra Z2 symmetry. Note that certain eigenvalues of simple and super hexagons 
for the degenerate Th bifurcation problem depend on the coefficients of quartic terms (see Table 13). 
These eigenvalues are now determined at quintic order. 

The stability results for the axial planforms are summarized in Table 16. The quintic truncation 
of the normal form is 



Note that the leading order a, (3 dependent terms are 0(2a — 1), where 2q! — 1 > 5, with 2q! — 1 = 5 
only in the case of (a, (3) = (3, 2). 

We assume that the following nondegeneracy conditions are satisfied: 



ii = Xzi + ziiailzil"^ + a2\z2\'^ + a2\z3\'^ + ttilzil"^ + aslzsl"^ + aelzel"^) 



+^i(/iiki|'* + f\2\z\^\z2^ + • • • + /seksPkeP + /eel^eP) 



(5.33) 



■\Z2Z-i{dxZxZ2Z-i -I- d2ZxZ2Z-i -I- d-i 

Z4Z5Z6 ~\~ '^4^4^5^6) 



+eiz-^ z^z^Zq + e2Zi z^ z^ z^ -\- ) . 



(ai + 02) 
(ai + 2a2) 
(ai - 02) 
2(ai - a2 f 



0, ±a25 ±a4, ias, i^e , 



(04 05), (04 ae), (as + a&) , 
0, ±(a4 as + ae). 



±(a4 - as), ±(a4 - ae), ±(as - ae) , 
(a4 - as)^ (a4 - ae)^ (as - a&f , 



(5.34) 




a /3 

0, ±^3, 

r if (a, /3) 7^ (3,2) 

I ±(^) if (a,/3) = (3,2) 



In this case we can draw a number of conclusions from Tables 13 and 16. 



1. The rhombs Rh/io always bifurcate unstably [13]. 
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Table 15: Eigenvalues for axial planforms associated with 12-dim. representations of Th + Z2. Also 
see Table 12. 



Axial Planform 



Eigenvalues 



Rhombs (Khho] 
1,0, 0,0,0) 



££l _|_ ££l (a) ££l _ ££l (a) £s£ dgl 

dxi dx2 ' dxi dx2 ' 9^3 ' dyl^ ' 

g (mult. 2), Ig (mult. 2), |g (mult. 2), (mult. 2) 



Simple Triangles 
y(l!,l!,l!, 0,0,0) 



9£l 
dy 



1^2M(a) M _ M (a) . 2), 

1 0|/2 ' dyi dy2 V /' 



g (mult. 6), 3M, (mult. 2) 



Anti- Hexagons 

:i, 1,1, -1,-1,-1) 



Same as Super Hexagons. 
See Table 12 



Super Triangles 
y{i,i,i,i,i,i) 



dg\ _ 




1 M 


, M _ 


J_ M _ 


L M 


M J 


lM , 


_ M _ 


M _ 


M _ 




dyi 


dy2 




9|/4 


9|/5 


9|/6 


9» ^ 


^ 9|/2 ^ 


9 1/3 


9 1/4 


9 1/5 


dye ' 


9a[ 




9a[ 




9P[ 


93[ 


9a[ 


93[ 




9a[ 


9a[ 


9a[ 


dxi 


8X2 


8X3 


1 9a;4 


9a;5 


9a;6 


dxi 


9a;2 


9a;3 




9a;5 


9a;6 ' 



9a;2 



9a;3 



(mult. 2), (mult. 2), //i, (mult. 2) 



9|/i 



9 1/2 



9£j_ 9£j_^ 2 

9|/4 dye , 



9|/i 
9 1/5 



9£| 

9 1/3 



93 



9 1/2 



9 1/3 



9 1/4 



9£| 

9 1/6 



9,J } , 



' dyi 



dg\ 
9 1/2 



9 1/3 



Anti- Triangles 
y{i,i,i,-i,-i,-i) 



Same as Super Triangles. 



Here the effect on Dg of the hidden symmetry t^^ (3-1) is included. 
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Table 16: Stability results for the hexagonal lattice bifurcation problem with F = F/^ + Z2, from 
Tables 12 and 15, and equations 5.6 and 5.33. For rolls and other rhombs, see Table 13. 



Axial 
Planform 


Signs of Nonzero Eigenvalues 


Branching Equation (x G R) 
and z G Fix(S) 




sgn(ai + 02), sgn(ai — 02), —sgn(ai — 02), 
sgn(a4 + uq — ai — 02), sgn(a4 + — ai — 02), 
sgn{a^ -\- — ai — 02) 


Aa; + (ai + a2)a;3 + 0{x^) = 0, 
and z = (s, a;, 0, 0, 0, 0) 


SiH 


sgn(ai + 202), sgn(ai — 02), 
sgn(a4 + + uq — ai — 202), —sgn{d2) 


Aa; + (ai + 2a2)a;2 + O(x^) = 0, 
and z = (s, a;, a;, 0, 0, 0) 


SiT 


sgn(ai + 202), sgn(ai — 02), 
sgn(a4 + + uq — ai — 202), sgn{d2) 


Aa; + (ai + 2a2)a;3 + O(x^) = 0, 
and z = (ia;, ia;, ia;, 0, 0, 0) 


Common 

to 
SuH„^^, 

SuT„^^, 
AT„,^ 


sgn(ai + 2a2 + (24 + (25 + ae), 
sgn(ai + 2a2 — 04 — as — ae), 
sgn(pi + 112) = sgn(ai - 02), 
sgn{jiiii2) = s^ra[2(ai — 02)^ — (a4 — as)^ 

-(04 - aef - (as - aef], 
and signs of remaining eigenvalues given below. 


Aa; + (ai + 2a2 + (24 + as + ae)^;"^ 
+ 0{x^) = 0, 
and z given below. 


SuH„^^ 


-sgn[(2a - (3)ei + (a + /3)e2], -sgn(d2 + ^3), 
-sgn[3id2 + (^4) - (2/3 - a)(ei - e2)a;2(«-3)] (») 


^7 ('7' 'T' 'T' 'T' 'T' 'T'l 

— ^ • c ^ • c ^ • c ^ • c ^ • c ^ • c y 


AH„^^ 


(-l)«+is5n[(2a - /3)ei + (a + /3)e2], -sgn{d2 - ds), 
-sgn[3(d2 - (^4) - (-l)"(2/3 - a)(ei - e2)a;2(«-3)] (») 


^7 1 '7' '7' '7' '7' '7' '7' 1 

^ I • c ^ • c ^ • c ^ • c ^ • c ^ • c y 


SuT„^^ 


(-l)«+is5n[(2a - /3)ei + (a + /3)e2], S5ra(rf2 - (^3), 
S5ra[3(rf2 + (^4) + (-l)"(2/3 - Q!)(ei - e2)a;2(«-3)] (») 


2 — ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ tZ/ ^ I' ^ 


AT„,^ 


-s£to[(2q! - /3)ei + (a + /3)e2], sgn(d2 + ds), 
sgn[3(d2 - (^4) + (2/3 - a)(ei - e2)a;2(«-3)] (») 


2 — ^ ^ tZ/ ^ ^ ^ ^ tZ/ ^ % ^ ^ tZ/ ^ ^ tZ/ ^ 



Note that the (a, /3)-dependent terms can be neglected here for all cases except (q!,/3) = (3,2). 
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2. If simple triangles and simple hexagons are neutrally stable at cubic order, then one and only 
one of the two branches is stable. The relative stability properties of these two solutions is 
determined at quintic order. 

3. It is possible for super hexagons, anti-hexagons, super triangles and anti-triangles to be 
unstable, even if they are all neutrally stable at cubic order. 

4. It is possible for all of the axial planforms to bifurcate supercritically, but none be stable. 

5. If any axial solution branch bifurcates sub critically, then rolls, super hexagons, super triangles, 
anti-hexagons and anti-triangles are all unstable. 

6. If rolls or super hexagons bifurcate sub critically, then all axial planforms are unstable at 
bifurcation. 

7. If simple hexagons and simple triangles are the only axial solution branches to bifurcate 
sub critically, then it is still possible that one, but not more, of the rhombs solutions is stable. 
Similarly, if rhombs Khhj^a,fj (j = 1,2, or 3) is the only axial solution branch to bifurcate 
sub critically, then it is possible for simple hexagons (or simple triangles) to be stable, or for 
one or more of the remaining rhombs solutions to be stable. 

8. If simple hexagons and the rhombs Rh/io bifurcate sub critically, then it is possible of one, 
but not more, of the other rhombs to be stable. However, if simple hexagons and one of the 
rhombs other than Rh/io bifurcate sub critically, then all axial solution branches are unstable. 

9. If two of the rhombs Khhj^a,fj (j = 1,2, 3) solution branches bifurcate sub critically, then it is 
possible that one, and only one, of the following solutions is stable: the remaining rhombs 
Khhj^a,fj, simple hexagons or simple triangles. However, if Rh/io and one of the other rhombs 
bifurcate sub critically, or if all three of the rhombs Khhj^a,fj solution branches are subcritical, 
then all of the axial solutions are unstable. 

10. If a is odd, then the only solution branches that can co-exist stably are simple hexagons (or 
simple triangles) and the rhombs Khhi^a,fj, ^^h2,a,fj, ^^h3,a,fj- Any combination of two of 
these states can bifurcate stably. It is also possible for all three types of rhombs to be stable 
simultaneously. However, if two or more of the rhombs are stable, then simple hexagons and 
simple triangles are unstable. If a is even, then it is also possible for super hexagons and 
anti-hexagons to be stable simultaneously, or for super triangles and anti-triangles to both be 
stable. 



6 Conclusions. 

We have investigated steady, spatially-periodic planforms which bifurcate from a spatially-uniform 
time-independent solution of _E(2)-equivariant and _E(2) -|- Z2-equivariant PDEs. We have done this 
within the framework of finite-dimensional equivariant normal forms of steady state bifurcation 
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problems. Our analysis applies on a center manifold associated with a doubly-periodic solution 
space of the PDEs. We considered separately the cases where the solutions are doubly-periodic on 
a square lattice and on a hexagonal lattice. 

In the case of generic bifurcation problems with De-j-T^-symmetry, a result of Ihrig and Gol- 
ubitsky [15] ensures that all of the axial solutions on the hexagonal lattice bifurcate unstably due 
to the presence of a quadratic term in the normal form. This means that not only rolls and simple 
hexagons are unstable at bifurcation, but also a countable set of rhombs are unstable. In order to 
capture stable solutions within a local bifurcation analysis, we considered a degenerate bifurcation 
problem in which the coefficient of the quadratic term is zero. We also considered the case in which 
an extra reflection symmetry ensures that no even terms are present in the normal form. 

For both lattices we determined the stability of the planforms which are guaranteed to bifurcate 
from the trivial solution by the equivariant branching lemma [12, 27]. In order to do this, we derived 
the normal form of each bifurcation problem. An order 2(a + (3) — 1 truncation of the normal form 
is required to completely determine the signs of the eigenvalues for all axial planforms periodic on 
a square lattice. In the case of the hexagonal lattice, an order (2a — 1) truncation is necessary. 
However, an important practical consideration is that much is already determined at cubic order. 

Previous studies focused on the "small box" limit for which the size of the periodic domain 
coincides with the wavelength of the instability; this leads to a bifurcation problem on for the 
square lattice and a bifurcation problem on C"^ for the hexagonal lattice. We have used and 

representations respectively for the symmetry groups associated with the square and hexagonal 
lattices. These apply when the periodicity of the lattice is much greater than the wavelength of 
the instability. This analysis extends the results of the earlier and C"^ bifurcation studies, both 
enlarging the number of planforms which are supported by the lattice and allowing for a wider 
class of disturbances in the stability analysis. For example, by considering all of the irreducible 
representations the stability of rolls, simple squares, simple hexagons, and simple triangles, to 
an infinite number of perturbations, can be determined. In particular, our bifurcation analysis 
provides a framework for addressing the relative stability of simple hexagons (or simple squares) 
and a countably-infinite set of rhombs. This is of particular interest in light of recent laboratory 
experiments on chemical Turing patterns in which a transition from a hexagonal pattern to a 
rhombic one is observed [14]. Our unfolding of the degenerate bifurcation problem on the hexagonal 
lattice provides a mathematical setting for investigating such a transition. 

There is no general bifurcation theoretic framework for computing the relative stability of 
squares and hexagons since no lattice supports them both. However, we are able to compute the 
stability of hexagons relative to rhombs, which are "almost square", i.e., which are composed of 
rectangles with aspect ratio that is close to 1 (see Table 9). For example, for the representation 
of De-i-T^ with {a,f3) = (4,3), the rhombs Rh/,3,4,3 are made up of rectangles with aspect ratio 
approximately 0.96; the angle between the wave vectors Ki and Ke in this case is about 92°. 

By not requiring the periodicity i of the lattice to coincide with the wavelength of the instability 
1/kc, we were able to investigate axial solution branches with periodicity 1/kc and simultaneously 
solution branches that have fundamental periodicity i ^ l/kc- We called the latter states su- 
per squares, super hexagons, super triangles, and anti-squares, anti-hexagons, anti-triangles. For 
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-E(2)-equivariant PDEs, there is an infinite family of these solution branches that is parameter- 
ized by an integer pair (a, (3). For the square lattice, the periodicity of these axial planforms is 
i = ^/a"^ -\- f3'^/kc, while for the hexagonal lattice i = \/ o? -\- f]"^ — af3/kc. The wavelength of the 
instability 1/kc determines the scale of the internal structure of these super and anti-states, while 
i determines their periodicity. By increasing a and (3 we obtain axial planforms that are periodic 
on larger and larger scales, all of which bifurcate from the trivial solution at A = 0. This is per- 
haps interesting in light of recent hydrodynamic experiments on quasi-patterns [10]. We emphasize, 
however, that the existence of a center manifold in the quasi-periodic case has not been established. 
Thus our analysis, which assumes that there is a center manifold, breaks down in the limit that i 
goes to infinity. 
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7 Appendix. 

In this appendix we give some of the details of our calculations of the isotropy subgroups that are 
presented in section 3. The computations are divided into two parts. In the first part, we compute 
the isotropy subgroups S of F which are translation free, i.e., for which S fl = {Id}, where Id 
is the identity element in F. In the second part, we consider the isotropy subgroups of F which are 
not translation free. 

The computation of the translation free subgroups of F is simplified by the following observation. 

Fet 

nH:F^H (7.1) 
be the projection of F into the holohedry H defined by 

(h,Q)^h. (7.2) 

(In the case that F = F/1-I-Z2, let IIjj : F D6-I-Z2.) The projection IIjj is a group homomorphism. 
If S C F is translation free then S is Iljj-isomorphic to a subgroup of H because ker(njj) C SflT^ = 
{Id}. Moreover, if Si and S2 are conjugate subgroups of F then njj(Si) and njj(S2) are conjugate 
subgroups of H. 

To find the translation free isotropy subgroups of F, we proceed as follows: 
1. We list all the subgroups G of H up to conjugacy. 
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2. For each subgroup G we then compute, up to conjugacy, all subgroups K of F that are 
isomorphic by IIjj to G. 



3. We determine which of the subgroups K are isotropy subgroups. 

For the isotropy subgroups S = Sz that are not translation free, we first show that certain 
components of z must be zero for it to possess a nontrivial translation symmetry. We then classify 
the isotropy subgroups associated with these points z. 

7.1 Square lattice case. 

In this section, we compute the isotropy subgroups for the eight-dimensional representations of 
F, = D4+T2. 



Translation free isotropy subgroups. 

Proposition 7.1 Up to conjugacy, the subgroups 0/ D4[i?^/2 7 ''"^^i] '^''^ 1; '^21-^^^], Z|[r2;j], Z2[T(i], 
Z4R^/2], Df[i2,,r,J, Bl[R^,Td] and B^R 

7r/27 '''xil- 

The proof follows from elementary group theory. 

Proposition 7.2 Up to conjugacy, the translation free subgroups ofT^ are 

1, Z^[i2.], Zf[r,J, Zi[Td], 

Zf[(r,,, (1/2,0))], Z4[i2^/2], mR^,T,,], Df[i2^,(r,,, (1/2,0))], 

Df[i2^,(r,,,(l/2,l/2))], Bi[R^,Td], D4[i2^/2, (r., , (1/2, 1/2))] , B^R^^^r,,] . 

The proof can be found in [8], where the translation free subgroups are called "shifted sub- 
groups". 

Proposition 7.3 The subgroups of Proposition 7.2 are all isotropy subgroups. 

Proof: To show that a subgroup is an isotropy subgroup, it is sufficient to find a point z G 
with symmetry S, i.e., to show that S = Sz = {7 G F^ : 7(z) = z} for some z G C"*. The conclusion 
of the proposition follows from 



D4 


= 1,1,1) 


D4 


= S(i, 1,-1,-1) 






= ^(1,2,1,2) 


Df 


= ^(1,2,-1,-2) 




Df 


= ^(1,2,2,1) 


Df 


_ 1 ^(1,2,2,-1) 
[ ^(1,2,-2,1) 


if a is odd 
if (3 is odd 


Z4 


= ^(1,1,3,3) 


zi 


= ^(8,28,8,-28) 




Zf 


= ^(8,28,-28,-8) 


Zf 


_ 1 ^(8,28,-28,8) 
[ ^(8,28,28,-8) 


if a is odd 
if (3 is odd 


z^ 


= ^(1,2,3,4) 


1 


= ^(8,28,38,48) 





The computations are simplified by the following observations: 
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1. The action of on is translation free. 

2. The action of D4 on permutes and/or conjugates the coordinates of z. 

Suppose that (h, 0)(z) = z for one of the z above. Because of 2, the only possible nontrivial action 
of on z G is to multiply some of the coordinates by —1. Hence 20 acts trivially on C*. By 1, 
this implies that 20 = in T^. Therefore we only need to consider (/i, 0) where = 0, ^i^i, ^£2 
or + and h G D4. I 



Non-translation-free isotropy subgroups. 

We now find the isotropy subgroups S C such that S fl 7^ {Id}. Our approach relies on first 
showing that if there is a nontrivial translation that acts trivially on some z G C"*, then at least 
two of the coordinates of z must be zero. This follows from the next lemma. 

Lemma 7.4 //" G acts trivially on three coordinates of z £ C*, then = G T^. 

Proof: It is sufficient to show that if acts trivially on Zj for j = 1, 2 and 3, then also acts 
trivially on z^. Because the action of on is translation free, it then follows that = G T^. 
We now show that 

K4 = aKi + 5K2 + cKs , (7.3) 

where a, b and c are integers. The conclusion of the lemma follows immediately from (7.3), since if 
Ki • 0, K2 • 0, K3 • G Z, then K4 -0 G Z. If we substitute the expressions for the Kj's given in 
Table 1 into (7.3), we find that 

(a + l)a - (b - c)/3 = 0, 

ia-l)a + ib + c)(3 = 0. (7.4) 
The following statements are equivalent: 

(i) There exist three integers a, b and c such that (7.4) is satisfied. 

(ii) There exist integers a, 5, c, r and s such that 

a + 1 = r(3 , a — 1 = sa , 

b — c = ra , b + c = —s(3 . 

(iii) There exist integers a, 5, c, r and s such that 

2a = r(3 + sa , 2b = ra — s(3 , 

2c = —ra — s(3 , 2 = r(3 — sa . 
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Statements (i) and (ii) are equivalent because a and (3 are relatively prime. Statements (ii) and 
(iii) are equivalent because the two systems of equations are equivalent. 

It follows from a well known result of number theory that if (a, (3) = 1, then there exist integers 
r' and s' such that 1 = r'(3 — s'a. Setting r = 2r' and s = 2s' proves (iii). I 

Lemma 7.4 ensures that the non-translation-free isotropy subgroups are, up to conjugacy, S = 
Sz7 where z satisfies one of the following conditions: 

1. z = 

2. z = (^1,0,0,0), ^1 G R, ^1 > 

3. Z = {zi,Z2, 0, 0), Zi,Z2 eB., Zi,Z2 > 

4. Z = (2:1, 0,2:3,0), Zi,Z3 eH, Zi,Z3> 

5. z = (zi, 0, 0, Z4), zi, Z4 £ a, zi, Z4 > 

Note that we have exploited the D4 symmetry to assume, without loss of generality, that zi 
whenever z 7^ 0. Moreover, the symmetry allows us to assume that the coordinates of z are 
real and nonnegative for the non-translation-free isotropy subgroups S = Sz. 

It is now a simple computation to determine, up to conjugacy, the non-translation-free isotropy 
subgroups listed in table 3. Presented below are the details of the computations for case (3). The 
others are similar. 

Let zi and Z2 be the only nonzero coordinates of z. The largest subgroup of acting trivially 
on zi and Z2 is 81^2 given in Table 3. We obtain 81^2 by solving 

Ki • = rai G Z and K2 • = ^2 G Z 

for G T^, where Ki and K2 are given in Table 1. We find that 

^ ( « f3_\ / -13 a \ 

The group 81^2 does not act trivially on any of the other coordinates because of Lemma 7.4 and 
the fact that 81^2 7^ {Id}. 

We used the action of T2 to assume that S = Sz where zi and Z2 are positive real numbers. 
Up to conjugacy, we obtain the following isotropy subgroups: 

S = Sz = Z^ + 81,2 if ^1 ^ Z2 

and 

S = Sz = Z4-|-8i^2 if ^1 = ^2 • 
Note that Z4 is the largest subgroup of D4 that preserves the set {±Ki, ±K2}. 



47 



7.2 Hexagonal lattice case: F = Th- 

In this section, we compute the isotropy subgroups for the twelve- dimensional representations of 
Th = D6+T2. 

Translation free isotropy subgroups. 

As for the square lattice case, if S C T/^ is translation free then S is 11 -isomorphic to a subgroup 
of De- Thus, to find all translation free isotropy subgroups of Th = De-j-T^, we first identify those 
subgroups of Th that are isomorphic by 11 to a subgroup of De- 

Proposition 7.5 Up to conjugacy, the subgroups of'De[R^^:j,Tx^] are 1, 7i2[R-k\, '^ii'^xi], Z2[t„], 
MRll^], MR^I^l Df[i2^,r,J, D3[i22/3,r,J, D-[i22/3,r„] anrf Deii^^/s, r^J- 

The proof follows from elementary group theory. 

Proposition 7.6 Up to conjugacy, the translation free subgroups of are the subgroups of De 
given in Proposition 7.5. 

The proof can be found in [8]. 

Proposition 7.7 The subgroups of Proposition 7.5 are all isotropy subgroups. 

Proof: The conclusion of the proposition follows from 





— ^(1,1,1,1,1,1) 


D3 


^(8,8,8,-8,-8,-8) 


Ds 




Ze 


= ^(1,1,1,2,2,2) 


Df 


= ^(1,2,3,3,2,1) 


Z3 


= ^(8,8,8,28,28,28) 




= ^(8,28,38,-8,-38,-2) 


Zf 


= ^(8,28,38,38,28,8) 




= ^(1,2,3,4,5,6) 


1 


= ^(8',28',38',48',58',68') 



The computations are similar to those for Proposition 7.3. In particular, if (h, 0)(z) = z for 
one of the z above, then 20 = in T^. Therefore we need only consider (/i, 0) where = 0, ^i^i, 
or + \l2 and /i G De. ■ 



Non-translation-free isotropy subgroups. 

We now find the isotropy subgroups S C T/^ such that S fl 7^ {Id}. We use the following lemma. 

Lemma 7.8 If Q £ acts trivially on four of the coordinates of z £ C^, then = G T^. 
When acts trivially on only three of the coordinates, then either (1) Q = £ T'^ , or (2) the three 
coordinates are zi,Z2,Z3, or they are 24,25, zq. 
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Proof: If acts trivially on two of the coordinates Zj, where j = 1, 2 or 3, then acts trivially 
on all three, because Ki + K2 + K3 = 0. Similarly, if acts trivially on two of the coordinates 
2:4, 2:5, zq then it necessarily acts trivially on all three because K4 + K5 + Ke = 0. 

We complete the proof by showing that if acts trivially on Zj for j = 1, 2 and 4, then also 
acts trivially on z^] thus it acts trivially on all z G C^. Hence = G since the action of 
on is translation free. 

Assume that acts trivially on zi,Z2,Z4^. It follows from the observation that if there exist 
integers a, h and c such that 

K5 = aKi + 5K2 + CK4 (7.5) 

then also acts trivially on z^. Substituting the expressions for the Kj's given in Table 2 into 
(7.5), we obtain the following conditions on a, 5, and c: 

(1 - 5 + c)q! + (a - c)/3 = 0, 

(a - 5 + c)a + (5 + l)/3 = 0. (7.6) 
The following statements are equivalent: 

(i) There exist three integers a, h and c such that (7.6) is satisfied. 

(ii) There exist integers a, 5, c, r and s such that 

1 — 5+c = r/3, a — c = —ra , 

a — h -\- c = s(3 , 5+1 = —sa . 



(iii) There exist integers a, 5, c, r and s such that 



a = r(p — a) — sa — 2 , b = —sa — 1 , 

c= r/3-sa-2, 3 = r(2/3 - a) - s(a + /3) . 



Statements (i) and (ii) are equivalent because a and (3 are relatively prime. Statements (ii) and 
(iii) are equivalent because the two systems of equations are equivalent. 

The proof proceeds by showing that (2/3 — a) and (a + (3) are relatively prime. To see this, 
suppose that d divides both 2(3 — a and a + (3, then d divides their sum 3(3. Since (3, a + (3) = 1, we 
know that d 3, so d must divide (3 as well as a + (3. Putting this all together we have d divides 
both a and (3, which are relatively prime so c? = 1. Hence (2(3 — a, a + /3) = 1 as claimed. It follows 
that there exist integers r' and s' such that 1 = r'(2(3 — a) — s'((3 + a). Setting r = 3r' and s = 3s' 
proves (iii). I 



Lemma 7.8 ensures that the non-translation-free isotropy subgroups are, up to conjugacy, S = 
Sz7 where z satisfies one of the following conditions: 

1. z = 
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2. 


z = 


(^1,0,0,0,0,0), 


zi e zi > 


3. 


z = 


(2:1, 2:2, 2^3, 0,0,0 


), zi,Z2 e R, 2:1,2:2 > 0, 2:3 7^ 


4. 


z = 


(^1,^2,0,0,0,0) 


, 2:1 , 2:2 G R, 2:1 , 2:2 > 


5. 


z = 


(^1,0,0,^4,0,0) 


, Zi, Z4 G R, ^1 , ^4 > 


6. 


z = 


(^1,0,0,0,^5,0) 


, 2:1 , 2:5 G R, 2:1 , 2:5 > 


7. 


z = 


(^1,0,0,0,0,^6) 


, 2:1 , 2:6 G R, 2:1 , 2:6 > 



It is now a straightforward computation to determine, up to conjugacy, the non-translation-free 
isotropy subgroups in table 4. We present the details of the computations for cases 3 and 4, only. 

Assume that 2:4 = 2:5 = 2:6 = and that zi and 2:2 are nonzero coordinates of z. The largest 
subgroup of acting trivially on 2:1 and 2:2 is 81^2,3 given in Table 4. We obtain 81^2,3 by solving 

Ki • = rai G Z and K2 • = ^2 G Z 

for G T^, where Ki and K2 are given in Table 2. It follows from Ki K2 K3 = that K3 • 
is also in Z. The group 81^2,3 does not act trivially on any of the other coordinates because of 
Lemma 7.8 and the fact that 81^2,3 7^ {Id}. 

Note that Zq is the largest subgroup of De that preserves the set {±Ki, ±K2, ±K3}. Hence, 
isotropy subgroups S = Sz for cases 3 and 4 are determined by considering, in turn, the action 
of (R^, 0), (i22^/3, 0), and (i2^/3, 0), G T^, on z = (n, r2, r3e2-'^, 0, 0, 0), where n, r2, r3 G R 
with ri,r2 > 0, > and G [0, 1). Here we have used the translation symmetry to assume 
that 2:1,2:2 are real and positive for some Sz in each conjugacy class. Moreover, we analyze cases 3 
and 4 simultaneously by allowing the possibility that 2:3 = in the following. 

1. (i?7r, 0)z = z. Substituting — Ki — K2 for K3, we have 

(i2^,0)z = (e-2-(Ki-0)^,,e-2-(K2-e)^^^g-2..c,g2..(Ki+K.).e^^^Q^Q^Q^_ 

Hence (i?7r, 0)z = z only if (/rf, 0) G 81^2,3 and either r3 = or (/j = 0, ^, i.e., 2:3 G R. 

2. (-R27r/3, 0)z = z only if ri = r2 = r3 = r. In this case, 

(i22V3, 0)z = (e2-.c.g-2..(Ki-e)^^ ^-2.^{K,.e)^^ ^2.^{K,+K,).e^^ (7_g^ 

and (-R27r/3, 0)^ = z provided we choose G such that K2 • G Z and Ki • — (/j G Z, 
which we can always do. 

3. Finally, we consider the equation (-R,r/3, 0)^ = z, where, again, ri = r2 = r3 = r and 

(i2^/3, 0)z = (e-2^''(Ki-e)^^ ^-2^^.^^-2MK2■@)^^ g2^.(Ki+K2)-e^^ 0^ 0^ ^7_9^ 

Thus G and (/j G [0,1) must satisfy the following equations: 

Ki-0gZ, K2-0±¥'GZ. (7.10) 
Hence (-R,r/3, 0)z = z for some G only if (/j = 0, ^. 



50 



From the above considerations we obtain, up to conjugacy, the following isotropy subgroups 
associated with the subspace {z G : = = zq = 0}: 

1. S = Sz = Si^2,37 where zi, Z2 and z^ are not all of the same norm, zi and Z2 are positive real 
numbers, and 2:3 G C \ R. 

2. S = Sz = Z2-(-Si^2,37 where zi, Z2 and z^ are not all of the same norm, zi and Z2 are positive 
real numbers, and z^ G R. (Note that this gives the isotropy associated with case 4 for which 

Z3 = 0.) 

3. S = Sz = Z3-i-Si^2,37 where zi = Z2 = z^ £ C \ Here we have used the observation that z 
is on the group orbit of (r, r, re^'^"^, 0, 0, 0), where r > is real and (f £ (0, 1), 7^ |. 

4. S = Sz = Z6-i-Si^2,3, where zi = Z2 = Z3 e R. 
7.3 Hexagonal lattice case: F = F/j + Z2. 

In this section we determine the additional isotropy subgroups that result from enlarging Th to 
Th + Z2. 

Translation free isotropy subgroups. 

Our approach is the same as for the previous two cases. We begin by finding, up to conjugacy, all 
subgroups of Th + Z2 that are isomorphic by Hj)g_|_22 to a subgroup of De + Z2. 

Proposition 7.9 Up to conjugacy, the subgroups of'De[R^^:j,Tx^] + Z2[k] are 

first class: 1, Zl[R^], Zf[r,J, Z-[r„], Z^[Rli^], Ze[R^,^], Df[i2^,r,J, T>^[RIi^,t,,,], 

second class: Groups of the form A + Z2 where A is one of the groups of the first class. 

third class: Z2[((i2., 0), k)], Z2[((r,, , 0), k)], Z2[((r„, 0), k)], Z6[((i2^/3, 0), k)], 
D2[((r., , 0), 0), iiR^, 0), K)], D2[((r,, , 0), k), ((i?^, 0), 0)], 
D2[((r., , 0), k), {{R^, 0), k)], BsiiiRlfs, «), 0), {{t,„0), k)], 
D3[((i2^/3, 0), 0), ((r„, 0), k)], B6[{{K/3, 0), 0), {{t,„0), k)], 
B6[{{K/3, 0), k), ((r.„ 0), 0)], and BeidK/a, 0), «), ((r.,, 0), k)]. 

The proof follows from elementary group theory. 

Proposition 7.10 Up to conjugacy, the translation free subgroups of Th + Z2 are the subgroups 
given in Proposition 7.9. 

The proof is similar to the proof of Proposition 7.5 given in [8]. 
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Proposition 7.11 Only the subgroups of first and third classes in Proposition 7.9 are isotropy 
subgroups. 



Proof: Because Fix(Z2) = {0}, we have Fix(S) = {0} for any subgroup in the second class. 
However, the origin has full symmetry F/^ + Z2, so no subgroup of the second class is an isotropy 
subgroup. 

For the groups in the first class, the conclusion of the proposition follows from 





= ^(1,1,1,1,1,1) 


Dg 


= ^(1+8,1+8,1+8,1-8,1-8,1-8) 


D3 


= ^(1+8,1+8,1+8,1+8,1+8,1+8) 


Ze 


= ^(1,1,1,2,2,2) 


Df 


= ^(1,2,3,3,2,1) 


Z3 


= ^(1+8,1+8,1+8,1 + 28,1 + 28,1 + 28) 


Z^ 


= S(l+8', 1 + 28', 1+38', 1-8', 1-38', 1-28') 


Zf 


= ^(1+8,1 + 28,1+38,1+38,1 + 28,1+8) 




= ^(1,2,3,4,5,6) 


1 


= S(l+8', 1 + 28', 1+38', 1+48', 1 + 58', 1+68') 



For the groups in the third class, the conclusion of the proposition follows from 

= Sr 



D6[((i2^/3,0),K),((r,,,0),K 
D6[((i2^/3,0),K),((r,,,0),0 
D6[((i2^/3,0),0),((r,,,0),K 
D3[((i2^/3,0),0),((r.„0),K 
D3[((i2^/3,0),0),((r„,0),K 
D2[((r.,,0),K),((i2^,0),K 
D2[((r.,,0),K),((i2^,0),0 
D2[((r.,,0),0),((i2^,0),K 

Z6[((i2^/3,0),K 

Z2[((r.,,0),K 
Z2[((r„,0),K 
Z2[k 



-'(«',«',«', — «', — ? 



-0 



(1,1,1, 



-1) 



(1+8,1+8,1+8,-1- 



-1-8) 



(1+8,1+8,1+8,-1+8,-1+8,-1+8) 
(8,28,38,-38,-28,-8) 



(1,2,3,-3,-2,-1) 

(8,28,38,38,28,8) 
(8,8,8,28,28,28) 

(1+8,1 + 28,1+38,-1-38,-1-28,-1-8) 
(1+8,1 + 28,1+38,-1+8,-1+38,1 + 28) 
(8,28,38,48,58,68) 



The computations are almost identical to those for the proof of Proposition 7.7. 



Non-translation-free isotropy subgroups. 

Lemma 7.8 still applies in the present context. We proceed, as in the case of F = F/^, by determining 
the isotropy of the following points in C^: 

1. z = 

2. z = (^1,0,0,0,0,0), ^1 G R, ^1 > 

3. Z = (zi,Z2,Z3, 0,0,0), Zi,Z2 e R, Zi,Z2> 0, Z3 ^ 0, 

4. z = (2:1, 2:2, 0,0, 0,0), zi,Z2 eB., zi,Z2 > 
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5. z = (zi, 0, 0, z^, 0, 0), zi, z^ £ a, zi, z^ > 

6. Z = (2:1, 0,0, 0,2:5,0), Zi,Z5 eH, Zi,Z5 > 

7. z = (2:1, 0,0, 0,0, 2:6), 2:1 , 2:6 G R, 2:1 , 2:6 > 

A straightforward computation determines the non-translation-free isotropy subgroups in table 
6. We present the details for case 4, only. 

Using the action of T^, we can assume that zi and Z2 are real. The group Z2[Rtt] + Z2[k] is the 
largest subgroup of De -|- Z2 that preserves the set {±Ki, ±K2}. The element i?^ acts trivially on 
z in this case, and k(z) = — z. We determine the generators of the new 81^2 by solving 



for G T^. We find that 81,2 is generated by ((Id, 2(a2-»fl+fl2) ^i - 2(«2-"«/h/32) ^2), and ((Id, 



K,.0 = -, 



and K2 • = ±- , 




2(a2_„^_^^2)*l 
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Figure 1: (a) Typical neutral stability curve, (b) Circle of critical wave vectors in k-space. 
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(a) 



Figure 2: (a) Critical circles for the square lattice when kc = 1, -\/3 and 5. (b) Critical circles for 
the hexagonal lattice when kc = 1 and \/l . 




Figure 3: Square lattice wave vectors Kj for the eight-dimensional representations of F^ in Table 1. 
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Figure 4: Hexagonal lattice wave vectors Kj for the twelve- dimensional representations of Th in 
Table 2. 




(d) (e) (f) 

Figure 5: Axial planforms associated with 8- dimensional representation of F^ with (a, (3) = (2,1) 
and xi,X2 G [—1, 1] (i.e., four copies of the fundamental domain are shown); (a) rolls, (b) rhombs 
(Rhsi^2,i)7 (c) rhombs (Rhs2, 2,1)7 (d) simple squares, (e) super squares, and (f) anti-squares. 
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(d) (e) (f) 



Figure 6: Axial planforms associated with 12- dimensional representation of Th with (a, (3) = (3,2), 
xi,X2 e [--^^-^h (a) rhombs (Rh/,1,3,2), (b) rhombs (Rh/,2,3,2), (c) rhombs (Rh/,3,3,2), (d) rolls, 
(e) simple hexagons (SiH"*"), and (f) super hexagons (SuH"*"). 
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(d) (e) (f) 

Figure 7: Axial planforms associated with 12- dimensional representation of Th + Z2 with (a, (3) = 
(3,2), xi^X2 G [— (a) rhombs (Rh/io), (b) super hexagons, (c) anti-hexagons, (d) simple 
triangles, (e) super triangles, and (f ) anti-triangles. (See Figure 6 for the additional axial planforms: 
rolls, simple hexagons, Rh/,1,3,2, Rh/12,3,2, and Khh3,3,2-) 
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Figure 8: Example of an hexagonal lattice bifurcation diagram for solutions in the six- dimensional 
subspace, z = (zi, Z2, z^, 0,0,0). Solid (dotted) lines indicate stable (unstable) solutions. The 
secondary solution branch has the form z = (si, a;2, a;2, 0, 0, 0), where xi,X2 G R. 
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Figure 9: Example of an hexagonal lattice bifurcation diagram for the twelve- dimensional repre- 
sentations of Th- Here < e ^ 1 in equation 5.29; see equation 5.30 for the other coefficients. 
Secondary bifurcation points are indicated by a solid circle; no secondary solution branches are 
shown. 




•. . ■. 
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